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ABSTRACT 

Aims. To investigate the spectroscopic properties of a selected optical photospheric spectra of core collapse supernovae (CCSNe). Special 
attention is devoted to traces of hydrogen at early phases. The impact on the physics and nature of their progenitors is emphasized. 

Methods. The CCSNe-sample spectra are analyzed with the parameterized supernova synthetic spectrum code “SYNOW” adopting some 
simplifying approximations. 

Results. The generated spectra are found to match the observed ones reasonably well, including a list of only 23 candidate ions. Guided 
by SN lb 19901, the observed trough near 6300A is attributed to Hff in almost all Type lb events, although in some objects it becomes too 
weak to be discernible, especially at later phases. Alternative line identifications are discussed. Differences in the way hydrogen manifests its 
presence within CCSNe are highlighted. In Type lb SNe, the Fla contrast velocity (i.e. line velocity minus the photospheric velocity) seems 
to increase with time at early epochs, reaching values as high as 8000 km s“' around 15 - 20 days after maximum and then remains almost 
constant. The derived photospheric velocities, indicate a lower velocity for Type 11 SNe 1987A and 1999em as compared to SN Ic 19941 
and SN 11b 1993J, while Type lb events display a somewhat larger variation. The scatter, around day 20, is measured to be ~5000 km s“*. 
Following two simple approaches, rough estimates of ejecta and hydrogen masses are given. A mass of hydrogen of approximately 0.02 Mq 
is obtained for SN 19901, while SNe 1983N and 2000F1 ejected ~0.008 Mq and ~0.08 Mq of hydrogen, respectively. SN 1993J has a higher 
hydrogen mass, ~ 0.7 Mg with a large uncertainty. A low mass and thin hydrogen layer with very high ejection velocities above the helium 
shell, is thus the most likely scenario for Type lb SNe. Some interesting and curious issues relating to oxygen lines suggest future investigations. 
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1. Introduction 

Stripped-envelope SNe, namely Type Ib (helium-rich) and 
Type Ic (helium-poor), being hydrogen deficient objects, 
are undoubtedly amongst the most mysterious SN classes. 
Recently efforts have started to understand the nature of these 
objects through studies of samples (Matheson et al. 2001; 
Branch et al. 2002). 

However the rarity of cases with well sampled observa¬ 
tions, photometry and spectra, hampers a more direct infer¬ 
ence of the physical situation behind the explosions and hence a 
clear view of the progenitor nature. Nevertheless, the discovery 
of metamorphosing events as SN 1987K and SN 1993J (recog¬ 
nized as SNe Ilb), that evolve from Type II to Type Ib-c as they 
age, together with the similarity of the environments in which 
they occur have linked Type Ib-c SNe to a core-collapse sce¬ 
nario in massive stars. At present, indeed, the most widely ac- 
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cepted models relate Type Ib-c SNe to both relatively low mass 
progenitors within the context of close binary system evolution 
(i.e. mass-loss as consequence of mass transfer) and massive 
stars that have undergone significant mass-loss due to a wind 
(i.e. Wolf-Rayet stars). So far observations have not discrimi¬ 
nated between the two scenarios, although recently, using HST 
data, a high spatial resolution search for the progenitor of the 
Type Ic SN 2004gt in a wide wavelength range from the far UV 
to the near IR has suggested that the event might result from an 
evolved Wolf-Rayet star, although the observations could not 
constrain models invoking less massive progenitors in binary 
systems (Gal-Yam et al. 2005). 

Photometrically, the lack of significant hydrogen in the 
outer layers of SNe Ib-c probably inhibits the most important 
characteristic of Type II SNe light curves, namely the plateau 
phase resulting from the hydrogen recombination wave. 

At late phases, the steeper decline rate, compared to the 
“^®Co to ^®Fe” decay slope, is indicative of significant y-ray 
escape as a result of the low mass ejecta in this class of ob- 
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jects (Clocchiatti & Wheeler 1997); there are rare exceptions 
where the late slope approaches the full trapping rate (e.g. SN 
Ib 1984L; Schlegel & Kirshner 1989). 

Spectroscopically, a clear separation scheme within the 
stripped-envelope SNe subclasses is still lacking. Part of the 
problem is the absence of meaningful statistics. A direct clas¬ 
sification was earlier proposed by Harkness et al. (1987) on the 
basis of He I strengths in the photospheric optical spectra of 
Type Ib SNe. Wheeler et al. (1994), however, have claimed the 
presence of He I 10830A in SN Ic 1990W, although He I lines 
were not noted in the optical region. The authors presented the 
idea of adopting, instead, the OI 7773A line as a distinguish¬ 
ing feature. The absorption seems stronger in Type Ic than in 
Ib SNe. Matheson et al. (2001) came to the same conclusion 
when analyzing a sample of Ib and Ic events. In addition, it 
has been argued that in Type Ib SNe He I lines 5876A and 
7065A gradually grow in strength with respect to the one at 
6678A (Matheson et al. 2001). SN 1998bw (GRB980425) is 
another example of classification as Type Ic, but where IR lines 
of He I have been clearly identified (Patat et al. 2001). 

Branch et al.(2002) have presented a relation between the 
velocity at the photosphere, measured using synthetic-spectrum 
analysis, and the time since maximum light for a sample of 
Type Ib SNe. The relation is well fitted by a power-law, in¬ 
dicating a fairly homogeneous behaviour. Interestingly, SN Ib 
19901 does not follow this trend (Elmhamdi et al. 2004). 

An important issue concerns features in the 
6000-6500A region of the early spectra in Ib-c SNe. 
Deng et al. (2000), when analyzing spectra of SN Ib 1999dn, 
have attributed the absorption feature seen at 6300A around 
maximum brightness to Ha which later on disappears or is 
overwhelmed by C II 6580A when Ha optical depth decreases 
as a consequence of the envelope expansion. The possibility 
of the presence of C II 6580A in Ib spectra has been ear¬ 
lier suggested by Harkness et al.(1987). Alternatively, the 
6300A absorption was identified to be due to Ne I 6402A in 
SN Ib 1991D (Benetti et al. 2002). Ne I lines and Ha had 
already been proposed to account for the deep absorption in 
SN Ib 1954A by Branch (1972). 

The presence, or not, of hydrogen and/or helium in Type 
Ib-c, with the possibility of quantifying the amount, is of great 
importance in identifying the progenitor stars that may give rise 
to these classes of objects. 

The papers by Matheson et al. (2001) and Branch et al. 
(2002) have provided an impetus towards an advanced under¬ 
standing of SNe Ib and Ic. The present work may be regarded 
as a continuation of those efforts. However in our comparative 
study of early spectra of Type Ib SNe we include representa¬ 
tives of all the various types of CCSNe, namely Type Ilb, Type 
Ic and Type II. This is established by means of synthetic spectra 
generated with the parametrized SN synthetic-spectrum code 
“SYNOW”. 

Our main goal was to understand the similarities and differ¬ 
ences among SNe Ib objects in the available sample, and to also 
compare to properties of the wider CCSNe family. The paper is 
organized as follows. First the analysis method and parameters 
are illustrated in Sect. 2. Data description is briefly given and 
the best fit synthetic spectra are presented and compared with 


the observed ones in Sect. 3. This is done separately for each in¬ 
dividual object among representatives of CCSNe classes. Sect. 
4 presents two methods to obtain spectroscopic mass estimates. 
The complete results will be discussed in detail and conclu¬ 
sions will be drawn in Section 5. 

2. Fitting procedure: SYNOW code 

For the purpose of our analysis we make use of the parameter¬ 
ized supernova synthetic spectrum code “SYNOW”. The code 
assumes spherical symmetry, homologous expansion: “v oc r” 
velocity-law and resonant scattering line formation above a 
sharp photosphere, emitting a blackbody continuum. SYNOW 
treats line formation and line blending (i.e. multiple scatter¬ 
ing) within the Sobolev approximation (Jeffery & Branch 1990; 
Fisher 2000; Branch 2001), and reads lines from the Kurucz 
42-million atomic line list (Kurucz 1993). It is an FTE code 
in only one respect: FTE excitation for the relative strengths 
of the lines of an individual ion, but it does not assume FTE 
ionization. 

SYNOW provides a number of free fitting parameters. The 
most important are: 1. (Tref), the optical depth of the strongest 
line, in the optical region, of the introduced ion. The line is 
called “the reference line”. The optical depths of the other lines 
of the same ion are thereafter determined assuming Boltzmann 
equilibrium. To decide which ion to introduce in the synthe¬ 
sis procedure we rely on the work by Hatano et al.(1999) as a 
starting point. The authors have presented the variation with 
temperature of FTE-Sobolev line optical depths of 45 indi¬ 
vidual candidate ions that might be encountered in supernova 
envelopes for six different compositions. 2. (Tht)' the underly¬ 
ing blackbody continuum temperature. We did not attach high 
physical importance to this parameter although the galactic ex¬ 
tinctions for each individual supernova, as reported by Schlegel 
et al. 1998, is taken into account. Independently of the present 
study, estimating the total reddening to the event, both that 
caused by foreground dust in the Milky Way and that caused by 
dust in the host galaxy, is a crucial point in supernovae study, 
especially when using them as cosmological probes. 3. (Vphot)- 
the velocity at the photosphere, estimated from the match with 
Fe II lines. Restrictions on the velocity interval within which 
an ion is introduced is possible. A maximum outer boundary 
velocity of line-forming-region of 5 x 10“* km s * is adopted 
in the present analysis. When assigning a minimum ion veloc¬ 
ity, “Vmin”, greater than “Vphot”, the ion is said to be “detached” 
from the photosphere, and consequently has a “non-zero” op¬ 
tical depth only starting at “Vmm”- In SYNOW the profile of 
a detached line has a flat-topped emission, and the absorption 
minimum is blueshifted by the detachment velocity. An unde¬ 
tached line has a rounded emission peak. A slightly detached 
line has a flat top but only over a small wavelength interval. 
However whenever one talks about detached lines one needs to 
keep in mind that we are using a somewhat unrealistic “Tv”, one 
that has a discontinuity in it. Real supernovae spectra probably 
do not have sharp discontinuities. 

The radial dependence of the line optical depths can be cho¬ 
sen to be either exponential with an e-folding velocity “v/’ as 
free parameter (i.e. t oc exp{—vIvg)), or a power-law with an 
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Table 1. The candidate reference lines of CCSNe 
phtospheric-phase spectra in the optical region (shown in 
increasing wavelength order). 


Ion 

Rest Wavelength (A) 

Ca II 

3934 

Nill 

4067 

Co II 

4161 

Mn II 

4205 

Cal 

4227 

CrII 

4242 

Sell 

4247 

CII 

4267 

Mg II 

4481 

Till 

4550 

Ball 

4554 

Fell 

5018 

Mg I 

5184 

He I 

5876 

Nal 

5890 

Si II 

6347 

Nel 

6402 

HI 

6563 

[Oil] 

7321 

OI 

7772 

Si I 

7944 

NI 

8680 

Cl 

9095 


index “n” (i.e. t oc v“"). In the present work a power-law pro¬ 
file is adopted with an index n = 8, although for some objects 
we will discuss the possibility of an exponential profile. 

When fitting our CCSNe sample spectra we have tested 
many combinations of fitting parameters. Only the best gener¬ 
ated synthetic spectra are displayed in the following and com¬ 
pared to the observed ones. We introduce the parameter “con¬ 
trast velocity”, defined as the line minimum velocity minus to 
thephotosphericone: “vcont(line) - Vmm(line)-Vphot”■ We show 
as well the behaviour of a similar parameter, defined instead as 
a ratio, i.e. = Vmin(line)IVphoi”- 

Although the computations are made under the purely res¬ 
onant scattering assumption, we find that the photospheric- 
optical spectra of CCSNe are fitted well usually requiring only 
23 candidate ions or fewer. Table 1 lists the candidate refer¬ 
ence lines sufficent to reproduce the observed features in opti¬ 
cal spectra of the CCSNe sample. It is important to note that 
since at early phases line formation takes place in high velocity 
layers causing severe line blending, it is better in some com¬ 
plicated cases to analyze line identifications in reverse chrono¬ 
logical order (i.e. starting with later phase spectra followed by 
earlier ones). All the analyzed spectra have been transformed 
to the rest frame of their host galaxies. 

3. Data description and Anaiysis 

Our selected CCSNe sample consists of 20 objects - 16 of them 
are Type Ib, 2 are Type Ic, 1 is Type Ilb and 1 is Type II 
SN. Some of the spectra presented in this work were gathered 


with the 60” Telescope and the “MMT” on Mount Hopkins. 
They are part of the supernova monitoring program of the 
Center for Astrophysics (PI R.P. Kirshner). Data for SN 1996aq 
are taken from the Padova-Asiago supernova database. Use is 
made as well of the “SUSPECT”' supernovae spectral archive. 
Descriptive data regarding the sample events are listed in Table 
2 (i.e. host galaxy, recession velocity and phases). For each in¬ 
dividual event, the table summarizes as well the most impor¬ 
tant fitting parameters. The last column indicates the number 
of ions we find responsible for determining the best fit spec¬ 
tra (details described in the following). Throughout the present 
work, however, we do not provide details concerning obser¬ 
vations of individual objects. We focus, instead, on synthetic 
spectra fitting procedures, constraints, problems and what we 
may learn about the CCSNe physical situation. 

3.1. Type Ib SNe: the sample 
* SN 19901: 

We start with SN 19901 since it can be considered as one of 
the better observed objects among Type Ib-c SNe in terms of 
quality and sampling. The event has been exploited both pho¬ 
tometrically and spectroscopically giving constraints on ejecta, 
oxygen and nickel masses and energy estimates (Elmhamdi et 
al. 2004). The supernova seems to follow a velocity trend dif¬ 
ferent from the pattern shown by the 10 SNe Ib sample of 
Branch et al. (2002). Here we investigate this peculiarity by 
means of synthetic spectra fits. 

The observed spectrum at maximum light, shown in Fig. 1, 
is compared to a synthetic spectrum (SSp) that has a velocity at 
the photosphere Vphot = 12000 km s“' and a blackbody contin¬ 
uum temperature Ti,b =14000 K. The SSp contains lines of He 
I, Fe II, Sc II, Mg II, O I, Ca II and Ha. The lines of Sc II, with 
t(Sc II)~ 0.4, have been introduced to help the fit of the absorp¬ 
tion feature around 6000A and redward the trough attributed 
to He I 5876A. Absorption troughs of P-Cygni He I lines at 
5876A, 6678A and 7056A are evident, although their relative 
strengths cannot be simultaneously fitted within the FTE ap¬ 
proach in the SSp. We will face this limitation each time we an¬ 
alyze and fit He I lines. A more precise analysis requires NFTE 
treatment as He I lines may be non-thermally excited by the 
decay products of ^®Ni and ^®Co (Fucy 1991). Apart from He I 
and Ha lines, the remaining lines have non-zero optical depths 
starting at the photospheric velocity (i.e. they are undetached). 
The SYNOW parameters required to account for the He I lines 
are T(He I)~ 2.9 and v™„(He I)= 14000 km s^'. 

The absorption minimum near 6250A is well fitted by Ha, 
with Vmin(Ha)= 16000 km s“* and assigned a moderate optical 
depth of 0.6 (Fig. 1). However, we tested other plausible alter¬ 
native identifications that were introduced in the literature to 
account for similar features seen in Type Ib SNe, namely Si 
II 6355A, Ne I 6402A and C II 6580A lines. Undetached Ne 
I lines were quite tempting in SN Ib 1991D, and worked as 
well as Ha having a contrast velocity of VcomiTta) =7000 km 
s“* (Benetti et al. 2002). For SN 19901, instead, Vconti Ha) is 
only 4000 km s-1 which means that undetached Ne I line is 

' http://bruford.nhn.ou.edu/ suspect/indexl.html 
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Fig. 1. The SN 19901 spectrum, near maximum, compared 
with the “SYNOW” SSp (thin line). Lines that are responsible 
for the most conspicuous features are also shown. The region 
around the 6250A trough is zoomed in the window. 


excluded because it is too blue to fit the observed feature. Si II 
6355A is obviously also too blue. C II 6580A remains then a 
plausible alternative for Ha, since its rest wavelength is slightly 
redder than Ha and therefore has its contrast velocity as a free 
parameter. In the window of Fig. 1, we show synthetic spectra 
for both cases (Ha and C II). C II, with v,„,„(CII)= 17000 km 
s ' and t(CII)~ 0.003, provides a fit as good as Ha. But an ac¬ 
ceptable C II 6580A line would require an additional velocity 
of about 820 km s * compared to Ha (~ 18A difference). C 
II 6580A has been proposed in early spectra of SN Ib 1999dn 
(Deng et al. 2000). Indeed, at maximum light and day 14, SSp 
fits of SN 1999dn included CII with v,o„,(CII)= 7000 and 1000 
km s ' respectively, while He I lines were formed starting at the 
photosphere (i.e.Vconf(He I)=0 km s ') (Deng et al. 2000). For 
SN 19901, the situation is almost similar; He I in the SSp shown 
in Fig. 1 has Vco«/(He I)=2000 km s"' while Vcon;(CII)=5000 
km s '. This means that the C II lines would need to be formed 
above the helium layer, which would be surprising, if not un¬ 
physical, in Type Ib SNe. 

SN 19901, with its particular velocity structure, presents 
good evidence for an Ha feature. However since it is not abso¬ 
lutely clearcut that Ha is always responsible for the absorption 
at about 6300A seen in Type Ib SNe, the presence of the Wfi 
Balmer line would be of great support. Unfortunately the opti¬ 
cal depth sufficient to fit the Ha trough is so small that the other 
Balmer features are too weak to be unambiguously detectable 
(T(Ha) is about 7 times greater than t(HjS)). 

Based on our investigation, the ions that generally might 
be encountered in shaping the 6000-6500A wavelength range 
are: Ha, Ne I, C II, Si II, Sc II, Ca I, He I, Fe II, Si I and Ba II. 



4000 5000 6000 7000 8000 9000 

Wavelength (A) 


Fig. 2. SN 19901 “SYNOW” fits of the observed 12day and 
21day spectra. Conspicuous line features are indicated 

Here, we propose the following criteria in view of the method¬ 
ology for deciding what lines should be adopted (especially for 
the 6300A feature): 

1/- when Fe II lines are very strong, they could produce a 
trough that might be sufficient to fit the 6300A feature. The 
depth of that Fe II feature is controlled by means of Fe II lines 
at 4924, 5018 and5169A. 

2/- undetached Sc II, Si I, Ca I and Ba II can be introduced 
especially to fit the slope at 6000-6300A wavelength range. 
They should not, however, introduce unwanted features in the 
rest of the spectrum. 

3/- undetached Ne I 6402A line is rejected once its feature is 
too blue to fit the 6300A trough or/and when the other Ne I 
lines clearly introduce various unwanted features. Similar rea¬ 
soning applies to the Si II 6355A line. 

4/- with its contrast velocity as a free parameter, C II 
6580A line could be a candidate for the 6300A trough; never¬ 
theless it is ruled out once it exceeds the He I contrast velocity. 

In Figure 2, the observed spectra at 12days (top panel) and 
21 days (bottom panel) are compared to the synthetic ones with 
(Vphot-^0000 km s“'; Thh-5500 K) and (vphot-9500 km s"'; 
Ti|(,=5400 K), respectively. Synthetic line profile features are 
labeled by the designation of the ion whose line gives rise to 
the feature. He I is still detached from the photosphere with 
Vco„,(He I)= 3000 km s^' at 12d and Vco„,(He I)= 2500 km s^' 
at 21days. The He I 7065A seems to increase in strength rela¬ 
tive to the He I lines at 5876A and 6678A, indicating the non- 
thermal excitation effects are changing but still existent. Two 
interesting points emerge from our SSp fit; on the one hand, 
even though the reference line Ca II 3933A, with T(Ca II)=120 
at both 12days and 21 days, produces a good match to the ob¬ 
served one, the Ca II infrared triplet (8542, 8662, 8498A) has 
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4000 5000 6000 7000 

Wavelength (A) 


Fig. 3. Synthetic “SYNOW” fit of the 21days observed spec¬ 
trum, with only lines of Sill (top), Nel (middle) and CII (bot¬ 
tom). 


a clear dehcit especially in the emission component of its P- 
Cygni profile. On the other hand, the observed O I 7773A is 
deep. An optical depth of 2 is in fact imposed to reproduce 
the deep trough in the SSp at 21days. The O I 7773A lines 
are believed to have relatively greater strength in Type Ic SNe 
compared to lb objects (Wheeler et al. 1994; Matheson et al. 
2001). This is presumably because for a “bare" C/O envelope 
of Type Ic, one would expect oxygen lines to be more promi¬ 
nent relative to Ib case where an intact helium, and possibly 
some hydrogen, could tend to dilute the C/O core. 

Matheson et al. (2001) have dehned a parameter called 
“Fractional Line Depth” of the line through absorption min¬ 
imum relative to the continuum flux. Mean values of 0.27 
(+0.11) and 0.38 (+0.091) were found to represent Type Ib and 
Type Ic respectively. For the case of SN 19901, we measure 
a value of about 0.45, indicating that the SN was peculiar in 
this respect. One explanation for this abnormal behaviour of 
SN 19901 might be related to the amount of oxygen. In fact 
Elmhamdi et al. (2004) have argued for a possible high oxygen 
mass (~0.7-l.35 Mq) relative to other Ib objects (~0.3 Mq in 
SNe 1984L, 1985Fand 1996N). 

As far as the 6000-6500A wavelength region is concerned, 
we hnd that Fe II, He I and Ha are sufficient to reproduce the 
overall shape in the SSp at both 12days and 21days. We have, 
however, tested the Ne I, Si II and C II possibilities as candi¬ 
dates for the 6250A trough. Panels in Fig. 3 are similar to the 
bottom of Fig. 2, but with only Si Il(top), Ne I(middle) and C 
Il(bottom) lines. All three ions are formed starting at the photo¬ 
sphere {Vcont- 0 km s '). Both Si II and Ne I lines are ruled out 
by means of criterion number 3 above. C II is as well rejected 
because of criterion 4. We therefore consider Ha to be the most 
likely explanation in SN 19901. 

* SN 1983N: 



4000 5000 6000 7000 8000 

Wavelength (A) 


Fig. 4. SSp fit of SN 1983N compared to the observed spectra 
at maximum light (upper panel) and at lOdays (lower panel). 
Conspicuous line features are shown. 

The observed spectra at maximum light and at ~10days 
are shown in Fig. 4, together with the corresponding synthetic 
spectra. The SSp at maximum brightness has Vphot - 11000 
km s“^ and a blackbody continuum temperature T*;, =8000 K, 
while the one at lOdays has Vphot - 7000 km s“' and T** =5000 
K. Ne I and C II are rejected because of the criteria 3 and 4, 
while undetached Si II helps in fitting the absorption blueward 
the 6250A trough. Ha with r ~0.34 accounts nicely for the 
6250A absorption at maximum. He I has t ~2.5 and Vconi- 0 
km s“', while SSp at lOdays has Vcont(He I)= 2000 km s“* and 
t'cof»(Ha)= 5000 km s'. Na I line is added in the lOdays SSp to 
help the ht of the emission component peaked around 5900A. 

* SN 1984L; 

Figure 5 compares the observed spectrum at ~8days with 
an SSp that has Vphot - 8000 km s“' and T^^, =8000 K. Ions 
that are responsible for the most conspicuous absorption fea¬ 
tures are indicated. Interestingly, He I lines with Vcom— 2000 km 
s * and T =2.5 , provide a good match at 5876A, 6678A and 
7056A. That may indicate that departures from LTE, due to 
non-thermal effects, are not severe for this Type Ib object. A 
similar result obtains for the 32day spectrum i.e. (Vco„,(He I)= 
2000 km s * and t =6.7; Top panel of Eig. 6). The correspond¬ 
ing SSp in Eig. 6 has Vphot - 5000 km s“' and T^,* =5600 K. 

The 32days He I profiles show rounded emission compo¬ 
nents that cannot be matched by a power-law SSp. As we have 
already mentioned, a detached line in a power-law assump¬ 
tion has a discontinuity in its optical depth (i.e. non-zero op¬ 
tical depth only above the VmmiUne)). When a profile retains 
a rounded P-Cygni emission component even if its VminiUne) 
is greater than Vphoi, this might indicate two components of 
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Fig. 5. SSp fit of SN 1984L compared to the observed spectrum 
at day 8. Conspicuous line features are shown. 


optical depth rather than only one above the photosphere. In 
such cases better fits could be obtained by having gradually 
decreasing optical depth below the detached velocity, instead 
of a discontinuity. In the present version of SYNOW, only 
an e-folding assumption of the optical depth allows a two- 
component treatment of a given line. We switch then to the 
exponential case adopting Ve=3000 (see Sect. 2). For He I lines 
and below Vmin- 7000 km s“* we introduce a second compo¬ 
nent with negative Vg (vg = -2000 km s“'), such that “r” is con¬ 
tinuous at the detachment velocity (i.e. at 7000 km s"'). In this 
case and for similar situations, the line should not be dehned 
as “detached”. Instead, it has a maximum value of “t” that is 
not at the photosphere as is normal for an undetached line. The 
bottom panel in Fig. 6 demonstrates that the two-component 
optical depth reflects more probably the real situation of the He 
I lines in SN 1984L. The He I line fits seem better relative to the 
power-law case, although blueward 5500A we obtain an infe¬ 
rior match. In the 8days SSp, Si II is too blue to account for the 
6250A. We checked the undetached Ne I 6402A line possibil¬ 
ity, assigning it an optical depth of 1.5. It gives a prohle broader 
than the observed one. In addition the coverage in wavelength 
of the spectrum hampers a check whether Ne I lines produce 
unwanted features longward of 7000A. Ne I lines would need 
to be non-thermally excited as the case for He I lines in Type 
Ib (Lucy 1991; Swartz et al. 1993). With an optical depth of 
1.5, Ne I would have a departure coefficient from LTE of about 
15 (Hatano et al. 1999). Additionally, the low non-thermal ef¬ 
fects seen in He I lines of SN 1984L (see above) may argue 
against Ne I identification. C II is also rejected since it will re¬ 
quire a Vco„,(Cll)= 8000 km s“' while at this phase Vcon/(He I)= 
2000 km s ' (criterion 4). Ha remains then a plausible candi¬ 
date. In Fig. 5, indeed, r(Ha)=0.65 and Veo„,(Ha)= 7000 km 



5000 6000 7000 

Wavelength (A) 


Fig. 6. SSp fit of SN 1984L compared to the observed spec¬ 
trum at day 32. The lower panel shows the e-folding optical 
depth possibility (see discussion). Conspicuous line features 
are shown. 



Fig. 7. SSp fit of SN 1988L compared to the observed spectrum 
at day 20. Conspicuous line features are shown. 


s * provide a good fit, while at 32days we use T(Ha)=0.2 and 
Vgg„,(Ha)= 7500 km s^' (Fig. 6). 


* SN 1988L: 
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Wavelength (A) 

Fig. 8. SSp fit of SN 1990B compared to the observed spec¬ 
trum around day 8. Lines that are responsible for the most 
conspicuous features are also shown. The region around the 
6300A trough is zoomed in the window. 



Fig. 9. SSp fit of SN 1990B compared to the observed spectrum 
around day 30. Lines of conspicuous features are shown. The 
region around the 5700A trough is zoomed in the window. Na I 
D and He I 5876A individual contributions are also displayed. 


Figure 7 compares the observed spectrum at ~20d with an 
SSp that has Vphot - 8000 km s ' and Tbh =9000 K. The 
observed spectrum has been smoothed with a box size of 5. 
Narrow emission lines due to H 11 regions are present in the 
spectrum. The object has been discussed and classified as a 
Type Ib SN (Filippenko 1988; Kidger 1988). 

The strong and broad P-Cygni profile at ~5550A is well 
matched in our SSp by a blend of Na I and He 1 lines. The 
He I 5876A accounts for the blue edge of the broad trough, 
having Vcon;(He 1)= 3000 km s ' and rjHe I)=l, while un¬ 
detached Na I has T(Nal)=3. Ba II lines together with Ca I 
lines have been introduced in the SSp and help the fit blue- 
ward ~6000A. Undetached Si II 6355A is too blue to account 
for the trough seen at ~6250A. We adopt Ha with an optical 
depth of 0.7 and Veo,!,(Ha)= 8000 km s '. Undetached C II 
6580A, t(CII)= 0.006, contributes together with He I blueward 
6500A (Fig. 7). SN 1988L may be regarded as an intermediate 
Ib/c object rather than a typical Type Ib. 

* SN 1990B: 

An extensive study of SN 1990B observations has been pre¬ 
sented by Clocchiatti et al. (2001). The authors pointed out the 
red character of the object. This fact is supported in the present 
analysis, see below, by means of the estimated continuum tem¬ 
peratures that seem much lower than values at similar phases 
for the other Type Ib-c sample objects. 

The spectrum at day 8 with the best fit SSp is presented in 
Figure 8. The SSp has Vphot - 8500 km s“' and Tbt =5400 K. 
Ions that are responsible for the most conspicuous features are 
indicated. The fit parameters have been modified many times in 


order to investigate line identification possibilities. The strong 
P-Cygni profile around 5700A is a blend of undetached Na 
I and detached He I 5876A (veon/(He I)= 1500 km s ' and 
T(He I)=0.55). The expelled helium, with Veo„,(He I)= 1500 
km s ', accounts nicely for the weak absorption near 6500A. 
Undetached Si II 6355A is considered to fit the weak absorp¬ 
tion at ~6200A. The Ha line, Vcont (Ha)= 4500 km s ' and 
T(Ha)=0.5, is introduced to account for the broad observed ab¬ 
sorption trough around 6300A. The window in Figure 8 dis¬ 
plays a zoom of the 6300A region, together with the C II fit 
possibility (t(CII)= 0.0025). The fit looks as good as Ha, how¬ 
ever it would need to be expelled with 3000 km s ' more than 
the He I lines (i.e. Vcon/(CII)= 4500 km s '), and hence is ruled 
out (criterion 4). 

The second spectrum, around day 30, is compared to an 
SSp that has Vphoi - 7000 km s“' and Tbb =4000 K (Fig. 9). 
The presence of He I 5876A is supported by two facts, first by 
the wide nature of the absorption trough that cannot be caused 
only by the undetached Na I D. Second, by the emergence of 
a “bump” around 5760A, in the transition slope redward of the 
deep absorption trough at ~5700A. A closer view of the “He 
I-i-Na I” region is provided in the window of Fig. 9. The win¬ 
dow demonstrates the contribution of both Na I D and He I 
5876A separately. The match is good and reproduces nicely the 
total feature seen in the observed profile. This tends therefore 
to confirm the presence of helium in the ejecta of this object, 
even though it does not exhibit a full and clear set of helium 
absorption lines. He I lines have VcoKr(He I)= 3000 km s ' and 
r(He I)=0.8. On the other hand, Fe II lines with T(FeII)=12, 
provide good fit in the blue part of the spectrum. In addition Fe 
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Days since explosion 

Fig. 10. The V-band light curve of SN 1990B compared to a 
sample of Ib-c “V-light curves” (See text for details). 

II contributes significantly in the absorption at ~6320A. The 
presence of Ha is therefore not needed since the Fe II contri¬ 
bution is sufficient. Nevertheless if Ha were present it would 
have a Vcon((Ha)= 4000 km s“'. 

Based on its spectroscopic properties, especially the weak 
He I lines, SN 1990B was re-classified as a Type Ic object 
(Clocchiatti et al. 2001). It was first classified as Type Ib SN 
(Kirshner et al. 1990). Clocchiatti et al. (2001) presented an 
extensive set of well-sampled photometry. In Figure 10, we 
show the absolute V-light curve of SN 1990B compared to 
other SNe Ib-c. The comparison highlights the similarity with 
Type Ilb 1993J, SN Ib 19901 and SN Ib 1991D. The maxi¬ 
mum brightness is comparable to those of SNe 1990 I, 1987M 
and 19941 and intermediate between the bright SNe 1999bw 
and 199ID and the faint SN 19831. Moreover the figure reveals 
a similar decline rate, from maximum to reach the exponen¬ 
tial decay, of SNe 1990B, 19901, 1993J, 199ID and 1998bw 
and indicates similar peak-to-tail contrast, whereas SNe 19831, 
1987M and 19941 display narrower peak widths and greater 
peak-to-tail contrast. The events with steeper decline rates (i.e. 
narrow widths) and greater peak-to-tail contrast are classified 
as Type Ic SNe, while SNe 19901 and 199ID are Type Ib 
events. This fact may point to a photometric behaviour closer 
to Type Ib rather than Type Ic. We regard, therefore, SN 1990B 
as an intermediate Type Ib/c event. 

* SN 1991ar: 

Figure 11 displays the available spectra of SN 1991ar at 
3 different phases, namely 28days, 48days and 60days. The 
first spectrum is compared to an SSp that has = 8000 
km s"' and Tth =4400 K. The broad P-Cygni profile with 
the minimum trough centered at ~5720A is well fitted by the 



Wavelength (A) 


Fig. 11. SSp fit of SN 1991ar compared to the observed spec¬ 
tra around 4 weeks (top panel), 48days (middle panel) and at 
60days (lower panel). Lines that are responsible for the most 
conspicuous features are reported. 

“He I-tNa I D” blend. Although the spectrum covers only the 
4000-7000A wavelength range, it shows distinctly the He I se¬ 
ries, indicating a Type Ib nature. The He I lines are undetached 
and have T(He I)=10. One particularity of this spectrum is the 
flat behaviour of the ~6070-6450A range. We checked differ¬ 
ent ion combinations, and the best fit to that zone is reached 
including lines of Fe II, Si II and Ha. This latter has Veo„,(Ha)= 
8000 km s ' and T(Ha)=0.32. 

For the 48day spectrum, the match is good with the SSp 
having Vphot = 7000 km s“' and Thb =4800 K (Fig. 11; middle 
panel). The He I lines are now slightly detached with VconAHe 
I)= 1000 km s ' and T(He I)=4.8 Lines of O I and Mg II 
both contribute to the observed features around 7600A and 
9030A (t(OI)= 3.5 and T(MgII)=5), while the Ca IIIR triplet 
is produced by T(CaII)=500. The weak feature labeled as be¬ 
ing due to Ha corresponds to VcontiHa)- 9000 km s ' and 
T(Ha)=0.4. 

The bottom panel in Figure 11 shows the 60day spec¬ 
trum compared to an SSp that has Vphor - 6000 km s“' and 
Thb =4600 K. The ~6080-6500A range loses the flat behaviour 
seen in the first spectra and a round emission profiles start to 
form. The match in that part of the spectrum with the SSp is 
acceptable. The profile is mainly due to lines of Fe II, Si II, C 
II, Ha and Ca 1. The He I fit corresponds to Vcon/(He I)= 1000 
km s“' and T(He I)=4.4, while Ha has Vco„f(Ha)= 8000 km s ' 
and r(Ha)=0.28. Note that the two weak troughs due to the Sc 
II lines at 5527 and 5661A appear near 5500A in all 3 spectra. 
Another notable feature is the deficiency, in flux, of the oh- 
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Wavelength (A) 

Fig. 12. SSp fit of SN 1991D compared to the observed spectra 
around 3 weeks. Lines that are responsible for the most con¬ 
spicuous features are reported. The top panel illustrates the Ne 
I possibility, while Ha case is displayed in the lower panel. 

served spectra, in the middle and bottom figures, relative to the 
synthetic ones especially redward of ~6500A. This is possibly 
because the supernova already enters a transition phase to the 
nebular epoch. 

Surely the evidence of Ha should not be taken as definite, 
although the corresponding weak absorptions in the SSp nicely 
matches the observed features. 

* SN 1991D: 

SN lb 199ID data has been presented and studied by 
Benetti et al. (2002). The object, with its narrow features, seems 
having lower velocities compared to other Type lb events at 
similar phases. When analyzing the ~10day spectrum by means 
of SYNOW code, the authors obtained an improved match in¬ 
troducing undetached Ne 1 lines with respect to the Ha SSp fit. 
Figure 12 illustrates this possibility for the 21 days spectrum. 
The displayed SSp has Vphot - 4600 km s“* and Ttb =7000 
K. The He I lines are evident through our SSp (Vcon/(He I)= 
1400 km s ' and T(He I)=1.8). In addition, the Na I D feature 
(T(NaI)=4) contributes to the He I 5876A broad P-Cygni pro¬ 
file. The presence of Na I lines is consistent with the good fit 
around 8100A. The most conspicuous line absorption features 
are indicated in the figure. In the top panel of Figure 12, un¬ 
detached Ne I lines are included in the SSp (T(NeI)=2). The 
match in the whole spectrum is quit good. In the bottom panel 
the Ha possibility is tested. Except for adding Ba II lines in or¬ 
der to help the fit at ~6100A, the other ion parameters are kept 
unchanged. The fit to the observed absorption near 6300A with 
Ha, VcontQio!)— 7400 km s“' and T(Ha)=0.46, is slightly better 
compared to the Ne I case (Top panel). However the SSp, in the 


bottom panel, does not account for the observed features near 
6630A and 6840A as does Nel lines in the SSp of the top panel. 
Ne I remains hence a strong candidate in this Type Ib object. 

It is worth noting here that at a similar phase, i.e. ~21days, 
SN Ib 19901 had a photospheric velocity of 9500 km s '. 

* SN 1991L: 

We study a spectrum dated about one month after max¬ 
imum (Fig. 13). The spectrum has been smoothed with a 
box size of 3. The classification of this event based only on 
the shown spectrum is very tentative. Indeed the spectrum 
shows a short wavelength part reminiscent of Type la spec¬ 
tra. However, nebular spectra of SN 1991L have been found to 
display emission features that are normally found in Type Ib/c 
events (Gomez & Lopez, 2002). The illustrated best fit SSp has 
Vphot - 5000 km s ' and T** =6000 K. The He I lines are 
not so obvious, although an undetached He I 5876A line, T(He 
I)=2, provides a good fit to the observed absorption trough at 
~5880A. Na I D could contribute to that feature, although it 
would provide an absorption component slightly redshifted. We 
regard SN 1991L as a transition Type Ib/c object. 

Similarly to SN 199ID, we tested the Ha and Ne I possi¬ 
bilities as a candidates for the feature near 6290A. Figure 13 
shows both cases. The included Ne I lines in the upper panel 
have T(He I)=l, while in the lower panel Ha has Vco«/(Ha)= 
8000 km s ' and T(Ha)=0.34. Some features due to Ne I lines 
appear similar to the observed ones, which are not accounted 
for the Ha case (lower panel; Fig. 13). Compared to the lower 
panel, Ne I lines also improve the fit of the broad emission 
component of the He I 5876A P-Cygni profile. Ne I remains 
therefore an alternative possibility to the Ha identification. 

* SN 1997dc: 

The best SSp fit, Vpi,ot - 5000 km s“' and T;,* =4200 K, 
is displayed in Figure 14 together with the observed spectrum 
(~4 weeks since maximum). On the one hand, the match with 
He I series, Veo„,(He I)= 2000 km s ' and T(He I)=2, is evi¬ 
dent, confirming the classification as a Type Ib event. On the 
other hand and in order to obtain an improved fit to the He I 
5876A P-Cygni profile, undetached Na ID with t= 5 is needed. 
The identification of Na I is consistent with the absorption pro¬ 
file seen near 8090A. Ca IIIR triplet is unusually strong and the 
match to the synthetic one is rather poor. Sc II with t= 2 is re¬ 
sponsible for absorption features blueward of the He I 5876A, 
and contributes as well, with Fe II, near 6160A. 

A zoomed view of the 5600-6600A region is shown in the 
window of Figure 14. The adopted SSp fits well the observed 
features. Various ion combinations have been tested, especially 
to explain the weak absorption at ~6290A. The best fit that 
would not introduce unwanted features in the rest of the spec¬ 
trum and as well be in agreement with our previous criteria 
(Sect. 3.1) is attributed to Ha having Vco«/(Ha)= 8000 km s“' 
and T(Ha)=0.2. 

* SN 1998dt: 

Spectra at 2 phases, ~8days and ~33days, are shown in 
Figure 15 and compared with the best fitting synthetic spectra. 
The SSp in the upper panel corresponds to vphot - 9000 km s“' 
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Fig. 13. SSpfitofSN 1991L compared to the observed spectra 
around 1 month. Lines that are responsible for the most con¬ 
spicuous features are reported. Top panel illustrates the Ne I 
possibility, while Her case is displayed in the lower panel. 

and Tift =5600 K, while for ~33days, the SSp has Vphot — 9000 
km s * and Tftft =5000 K (lower panel). The SSp nicely match 
the features in the observed spectra. The narrow emission near 
6565A is due to Ha from H II region. 

The reference Ca II line (i.e. at 3933A) in the upper panel 
has an optical depth of 500. The corresponding Ca II H&K 
and IR Ca II triplet (8542, 8662, 8498A) SSp both fit well the 
observed broad P-Cygni profiles. For the 33day spectrum how¬ 
ever, the Ca II infrared triplet (t=500) presents a deheit with 
respect to the observed prohle even though for the absorption 
part the ht is still acceptable. The O I 1113A with r=0.5 ac¬ 
counts for most of the observed feature at both phases (top and 
lower panels). 

The He I lines are clearly evident for both phases, with the 
main difference that while at ~8days He I lines are detached 
(v'con((He I)= 2000 km s * and T(He I)=4), they have non-zero 
optical depths starting at the photosphere for the ~33day spec¬ 
trum (vconAHe I)= 0 km s“' and T(He I)=10). Note here that 
for both phases He I 5876A is sufficient to fit the observed P- 
Cygni prohle. There is no need to include Na ID line. The slope 
redward of the He I 5876A emission component is nicely ac¬ 
counted for by lines of Fe II, Ca I and Ha at ~8days, while 
for the 33day spectrum we introduce, in addition, lines of Sc II 
(t=2). The presence of Sc II is supported by the double absorp¬ 
tion features near 5390 and 5500A. 

The best SSp ht of the observed 8day spectrum includes 
only 6 elements, namely Fe II, He I, Ca I, Ca II, O I and H I. 
Ha, with Vcont= 8000 km s ' and t= 0.3, accounts for the weak 
absorption near 6200A. A similar but weaker absorption is seen 
in the 33day spectrum as well. We ht it with Ha that has Vconi= 



Fig. 14. SSp ht of SN 1997dc compared to the observed 
spectrum at about 4 weeks. Lines of conspicuous features are 
shown. The region around the 6300A weak trough is zoomed 
in the window (see text). 

8000 km s ' and t—0.2. It is clear that even though the ht is 
convincing, the presence of Ha is not dehnite. This limitation 
occurs whenever we deal with absorption features that are not 
so deep with respect to the continuum. 

* SN 1998T; 

Spectra of SN 1998T at two different phases, ~ 21 days and 
42days, are displayed in Figure 16 and compared to synthetic 
spectra having Vpi,ot = 6000 km s“' and Tftft =5400 K (up¬ 
per panel), and Vphot = 5800 km s“' and Tftft =5200 K (lower 
panel). The spectra are highly contaminated by weak emission 
features from an H II region. The absorption troughs caused by 
He I series are prominent, indicating a familiar Type Ib appear¬ 
ance. The He I line best ht at 21 days corresponds to Veo„,(He 
I)= 3000 km s“' and T(He I)=L15, while at 42days Veo„f(He 
I)= 2200 km s ' and r(He I)=0.75. The optical depth of the 
undetached Ca II is 120 (upper panel) and 140 (lower panel). 
The corresponding features, i.e of Ca II, ht well the observed 
ones although an excess in the emission P-Cygni components is 
visible, as is the case in most advanced photospheric spectra. C 
I, with t= 0.3, helps ht the absorption feature near 8930A. Fe II 
lines are unusually weak for SN 1998T. Indeed, optical depths 
of 1 and 5 are respectively adopted for our best synthetic hts at 
~21days and ~42days. 

We checked different combinations for the 
6000-6500A wavelength range, testing the possible can¬ 
didates and keeping in mind our departure criteria (Sect. 3.1). 
For the 21 day spectrum, the best ht includes Si I(t= 0.03), 
Si II(t= 0.6), Fe II and Ha. This latter has Vcont (Ha)= 4500 
km s"' and T(Ha)=0.34. When testing Ne I lines one obtains 
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Fig. 15. SSp fit of SN 1998dt compared to the observed spectra 
at 8days (upper panel) and at 33days (lower panel). Lines of 
conspicuous features are indicated. 



4000 5000 6000 7000 8000 9000 

Wavelength (A) 


Fig. 17. SSp fit of SN 1999di compared to the observed spectra 
at 21days (upper panel) and at 45days (lower panel). Lines of 
conspicuous features are indicated. 
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Wavelength (A) 

Fig. 16. SSp ht of SN 1998T compared to the observed spectra 
at 3 weeks (upper panel) and at 42days (lower panel). Lines of 
conspicuous features are indicated. 

an improved fit blueward of the He I 7065A feature, but then 
unwanted profiles in the 6000-6500A region would be intro¬ 
duced. At low temperatures, an Si I identihcation is plausible 
in Type Ib objects within LTE assumption with low optical 
depths (Hatano et al. 1999). Si II however, is expected (in LTE) 


to have an optical depth of about 10 at temperatures similar 
to that for the 21 day spectrum (helium-rich composition; 
Hatano et al. 1999). This would mean that we might have 
mis-identihed the absorption feature near 6250A, or we have a 
departure from LTE in the Si II line. Note that for Type II SNe, 
and at similar temperatures. Si II has an optical depth as low 
as 0.1 (hydrogen-rich composition; Hatano et al. 1999). The 
remaining possibility is that the 6250A is due to Ha rather than 
Si 11. In this case hydrogen would have Vconf(Ha)= 9000 km 
s ', but then the absorption near 6340A remains unaccounted 
for. To fit this latter with C II, we need then to impose a 
velocity of about Veo„f=6000 km s“*, which means that carbon 
is expelled 3000 km s * greater than helium, contradicting our 
criteria (Sect. 3.1). The most acceptable combination is then 
the one illustrated in the top panel of Eigure 16. 

The situation for the 42day spectrum is more difficult as 
the spectrum is noisier. Our best ht includes Si I(t= 0.025), Si 
II(t= 0.4), Ee II and C II (t=2.3x 10-4). We prefer C II, having 
in this case similar Vcont as He I, rather than adopting Ha. Eirst, 
because hydrogen would need then to be expelled at velocities 
lower than helium, and second because we remain consistent 
with our criteria. We propose that it was Ha at 21 days which 
was then overwhelmed by C II later as the photosphere recedes 
because of the envelope expansion. 

* SN 1999di: 

In Eigure 17, the observed spectra of SN 1999di at ~21days 
and 45days are compared to the computed synthetic spectra. 
The spectra are synthesized adopting Vphot - 7000 km s“* 
and Tii, =4800 K (upper panel) and Vphot = 6000 km s“' and 
Thh =5200 K (lower panel). The match is good in the overall 
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Fig. 18. SSp fit of SN 1999dn compared to the observed spec¬ 
tra at lOdays (upper panel) and at 38days (lower panel). Lines 
of conspicuous features are indicated. The region around the 
6300A weak trough is zoomed in the window (see text). 


spectral shape. The fit to the He I lines is obvious. The pa¬ 
rameters used are Vco„,(He I)= 1500 km s ' and T(He I)=12.6 
for the 21 day spectrum, while at 45days the He I optical depth 
diminishes to T(He I)=10, keeping the same contrast velocity 
(i.e. 1500 km s '). For the 21day spectrum the optical He I 
lines (i.e. at 5876, 6678 and 7065A) are simultaneously well 
accounted for by their corresponding absorption troughs in the 
SSp. Even the He I 7281A line accounts for an observed ab¬ 
sorption. At day 45, however, the SN displays a shallower He 
I 6678A compared to the He I lines at 5876A and 7065A. At 
both phases, Na ID 5893A is blended with He 15876A. The ab¬ 
sorption near 8030A is also accounted for by lines of Na I. The 
features blueward of the broad “He I-rNa ID” P-Cygni prohle 
are attributed to Sc II lines (t= 2 for both epochs). Sc II also 
helps the ht redward of the emission component of the blended 
“He IH-Na ID” P-Cygni prohle. 

The O I 7773A line, with t= 1, accounts for most of the 
observed feature. Note however that the line of Mg II at 
7890A contributes in some cases to the red edge of the O I 
7773A line. In the present case this contribution is not so rel¬ 
evant. In the upper panel, the Ca II IR triplet, with an opti¬ 
cal depth of 700, hts very well the observed prohle. At day 
45, the ht is under-estimated especially for the emission part, 
which may indicate a transition to the nebular phase. The re¬ 
produced Ca II H&K line is also strong in both synthetic spec¬ 
tra. The issue of the 6300A feature is particular in this object, 
because it is very deep and different from features seen in “nor¬ 
mal” Type Ib SNe. We have tested various identihcations for 
the deep feature. The best ht that would be consistent with our 
criteria attributes the trough to Ha with Vcon/(Ha)= 5000 km 


s ' and T(Ha)=3 (upper panel), and Vco„f(Ha)= 6000 km s ' 
and r(Ha)=1.6. With these parameters and in addition to the 
almost perfect Ha ht, the SSp accounts as well for the notch^ 
seen near 4665A, assigning it to Wfi. The Ha trough is rem¬ 
iniscent of what is seen in Type II and Ilb objects, with the 
difference that in these latter types Ha displays a complete P- 
Cygni prohle (i.e. with both conspicuous emission and absorp¬ 
tion components). This is understandable within the context of 
“detachment” concept. Indeed, Branch et al.(2002) have shown 
that an Ha P-Cygni prohle loses its obvious emission compo¬ 
nent when it is highly detached. If correct, we would then ex¬ 
pect an increasing value of Vcon/(Ha) as we go from Type II to 
lib to Ib SNe. This point will be emphasized and discussed at 
the end of the paper. 

4 SN1999dn: 

Two spectra of SN 1999dn, observed at 10 and 38 days 
after maximum light, are compared in Figure 18 to synthetic 
spectra that have Vphot - 7000 km s“' and Tbh =5800 K (up¬ 
per panel) and Vphot - 6000 km s“' and T^i =5400 K (lower 
panel). As shown in both panels, the prominent lines of He I 
5876, 6678 and 7065A are clearly accounted for, indicating a 
typical Type Ib nature. The ht to the He I reference line corre¬ 
sponds to Vcon/(He I)= 2000 km s ' and T(He I)=1.9 at day 
10, and Vcon/(He I)= 1000 km s“' and r(He I)=14.5 at day 
38. A spectral analysis of SN 1999dn has been also presented 
by Deng et al. (2000). The authors discussed line identihca¬ 
tions for 3 photospheric spectra (at -10, 0 and 14 days from 
maximum). The observed trough at ~6250A was hrst blended 
with He I 6678A for the -lOday spectrum, while around max¬ 
imum light the two troughs become distinctly isolated, giving 
rise to a double absorption prohle (Fig. 2; Deng et al. 2000), 
very similar to what is seen in SN 19901 (Fig. 1), although 
the two events evolve afterwards in different ways, with the 
~6250A feature being less deep in SN 1999dnthan in SN 19901 
later on. The authors found it difficult to attribute the minimum 
near 6250A to Si II6355A. They argued that it was Ha hrst, be¬ 
fore maximum, that becomes blended and overwhelmed by C 
II line in later spectra. For the two early spectra, i.e. at -lOdays 
and maximum, C II 6580A provides a ht as good as Ha, how¬ 
ever C II 6580A is assigned a minimum velocity much higher 
than the one attributed to He I. Highly detached C II lines are 
surely improbable in this class of event. Even in their late spec¬ 
trum (i.e around 14days), the C II has a value of Vcom higher 
than that of He I. 

We analyzed various possibilities among line candidates for 
our 10 and 38 day spectra. The wavelength range of interest, 
~ 5800-6600A, is zoomed in the window of Figure 18. We hnd 
that the best hts which would reproduce the weak absorption 
features and not contradict the ht criteria, include lines of He I, 
Fe II, Ca I and Ha. This latter accounts nicely for the absorp¬ 
tion near 6260A with Vco«/(Ha)= 8000 km s“' and r(Ha)=0.24 
(at lOdays), and Vco«/(Ha)= 8000 km s“' and T(Ha)=0.34 (at 
38days). 


^ The word ’’notch” referes to a narrow and weak ahsorbtion fea¬ 
ture. 
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Fig. 19. SSp fit of SN 1999ex compared to the observed spec¬ 
trum at 4days. The observed IR spectrum is combined with the 
optical one in order to allow the He 110830A identihcation (see 
discussion). The Lines of conspicuous features are reported. 

* SN 1999ex: 

A good set of early data for SN 1999ex, spectra and pho¬ 
tometry, have been presented by Hamuy et al. (2002). The qual¬ 
ity of the photometry, starting well before maximum light, of¬ 
fers a unique possibility of looking at the early behaviour of 
light curves. Indeed, the early “dip" seen in the "U" and "B" 
light curves is interpreted as being due to the shock breakout, 
supporting the present belief that Type Ib-c SNe are the out¬ 
come of core collapse in massive stars rather than thermonu¬ 
clear disruption of white dwarfs (Hamuy et al. 2002; Stritzinger 
et al. 2002). Because of weak optical He I lines, the object was 
classihed as an intermediate case between Ib and Ic SNe. The 
evident trough at ~6250A was attributed to Si II6355A (Fig. 4; 
Hamuy et al. 2002). The authors also presented three infra-red 
spectra. 

We analyze, by means of synthetic spectra, line identifi¬ 
cations in this interesting object. The two observed spectra, 
around 4 and 13 days, are compared with our best ht SSp that 
have Vphot - 10000 km s“' and Thb =5800 K (Fig. 19) and 
Vpb„, = 7000 km s-' and Ttb =5600 K (Fig. 22). In Figure 19, 
we combined the IR spectrum with the optical one in order to 
check the consistency of He I identihcation. In fact, at 4days, 
assigning the following parameters: VconAHe I)= 1000 km s~' 
and r(He I)=2.35, we obtain a good match to the observed He 
I prohles. On the one hand, a strong support for the presence 
of He I lines comes from the good ht to the IR He I 10830A. 
Furthermore, in hgure 20 the ht is extended beyond 1 micron. 
The good match to other He I-IR lines, adopting equal v,„,„(He 
I) and r(He I) as in the optical part of the spectrum, is clear and 
provides unambiguous evidence for the presence of helium in 



Fig. 20. The 4day IR spectrum of SN 1999ex compared with 
synthetic spectra that contain only lines of He I. The prominent 
features attributed to He I lines are indicated. 
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Fig. 21. The 4day spectrum of SN 1999ex is compared with 
synthetic spectra that contain only lines of Ne I (top panel) and 
Si II (bottom panel). 

this event. The OI 7773A line with t=0.5, on the other hand, is 
found to be not so strong and deep as is the case in most Type Ic 
events (Matheson et al. 2001). These two facts point to a Type 
Ib nature, more than a Type Ic class. 

In the lower panel of Figure 19 we test the exponential case 
(see details in the “SN 1984L” part; Sect 3.1). We assign two 
components to the He I lines: one above Vmin- 11000 km s ' 
with Ve =3000 km s“* and one below Vmin- 11000 km s“* with 
negative Vg {Vg - -2000 km s“*), such that t is continuous 
at the detachment velocity (i.e. at 11000 km s“*). The ht is 























14 


Elmhamdi et al.: CCSNe spectra 



Fig. 22. SSp fit of SN 1999ex compared to the observed spec¬ 
trum at day 13. Lines of conspicuous features are shown. 

slightly improved for He I lines at 5876A and 10830A and also 
for the IR triplet Ca II. 

As far as the trough around 6250A is concerned, we 
checked the Si II identihcation attributed by Hamuy et 
al.(2002). We test as well the Ne I possibility. The closer 
view in Figure 21 illustrates the two possibilities. In the upper 
panel undetached Ne I lines can ht the feature, however they 
would introduce unwanted features in the rest of the spectrum. 
Furthermore, because of the depth of the observed trough, a 
good match with Ne 16402A would require an optical depth of 
4, indicating a departure factor from “LTE” of ~40 (Hatano et 
al. 1999), highly improbable although non-thermal excitation 
of Ne I needs to be more thoroughly investigated (i.e. through 
NLTE and hydrodynamic SSp codes). The Si II possibility is 
shown in the bottom of Eigure 21. The undetached Si II6355A, 
with T =2.5, is rather blue to account for the feature. We note 
here that if one adopts C II6580A, then one needs to assign it a 
very high velocity, about 8000 km s ' higher than the one of He 
I. The most likely identihcation therefore remains Her. In fact 
the best ht in Eigure 19 is achieved using Vconf(Ha)= 8000 km 
s * and t(Hq')= 1.45, while for day 13 Ha has Vconf(Ha)= 5000 
km s * and r(Ha)=1.6 (Eig. 22). The match at both phases is 
quite good. The other lines responsible for shaping the spectra 
are labeled in the hgures. In addition to the convincing match 
to lines of He I, IR-Ca II, O I, Ee II we note the presence of 
a weak absorption near 4680A that is well accounted for by 
Hj6 in our 13days SSp (Fig. 22). This can be taken as support 
for the presence of hydrogen expelled at high velocities in this 
supernova. 

It would have been of great interest to look at later spectra 
of this object, in order to see how the trough assigned to Ha 
would evolve with time. 
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Fig. 23. SSp fit of SN 2000H compared to the observed spec¬ 
trum at maximum light. Lines of conspicuous features are in¬ 
dicated. The top panel shows the “Ha-i-Ne I” possibility for the 
broad feature near 6300A. A corresponding zoomed view is 
displayed in the window. Similarly, the bottom panel shows the 
“Ha-tSi 11” case. 


* SN 2000H: 

The Type Ib SN 2000H is considered one of the more inter¬ 
esting Ib-c objects, especially with a strong and deep trough 
near 6300A (Benetti et al. 2000; Branch et al. 2002). Here 
we analyze three spectra at different epochs, namely at max¬ 
imum light, 19days and 30days. Figure 23 compares the ob¬ 
served spectrum around maximum with the resulting best ht 
SSp which has Vphot = 11000 km s * and Thb =10000 K. 
The overall match is quite good. Undetached He I lines, with 
T =2, provide a good ht with the He I 5876A, while features 
assigned to He I lines 6678A and 7065A are very weak. The 
absorption trough near 6280A is exceptionally broad and can¬ 
not be accounted for only by Ha. We identify two possibili¬ 
ties for which we obtain a broad absorption in agreement with 
the observed feature. The upper panel in Figure 23 shows the 
“Ha-i-Ne I” combination, while in the lower panel the “Ha-tSi 
11” combination is illustrated. In both panels a closer view of 
the 6280A region is displayed in the window. In the hrst case, 
undetached Nel with r =1.5, a contribution blueward of He 
I 7065A improves the ht with the observed feature, however 
the He I 5876A emission component is under-estimated. Note 
that the narrow absorption in the emission peak is attributed to 
Na ID interstellar line originating in the parent galaxy. In the 
lower panel, undetached Si II (r =3) improves the match with 
the observed absorption trough without altering the resulting 
SSp at the emission part of the He 15876A P-Cygni prohle. As 
a result we believe that the“Ha-i-Si 11” combination is the most 
probable. 
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Fig. 24. SSp fit of SN 2000H compared to the observed spec¬ 
trum at 19days. The lower panel shows the e-folding optical 
depth possibility (see discussion). Conspicuous line features 
are shown. 



Wavelength (A) 


Fig. 25. SSp fit of SN 2000H compared to the observed spec¬ 
trum at 30days. The lower panel shows the e-folding optical 
depth possibility. Conspicuous line features are indicated. 


Figures 24 and 25 compare spectra at 19 and 30 days with 
their corresponding best fit SSp. These have Vphot = 6000 km 
s“* and Tbb =4600 K (Fig. 24) and Vphot — 5000 km s“' and 
Tbb =6000 K (Fig. 25). The two spectra are similar in hav¬ 
ing naiTower features than the maximum light spectrum, re¬ 
sulting from observing at smaller radii where expulsion veloc¬ 
ities are lower. Lines of Ca II, [O II], O I, Mg II are seen to 
develop. Introducing Sc II lines helps to form a feature blue- 
ward of the strong He I 5876A P-Cygni profile. A distinct fea¬ 
ture appears redward of the IR Ca II profile, more evident in 
the 30day spectrum, and is accounted for by C I (r =0.2 at 
19days and t =0.8 at 30days). The noticeable change com¬ 
pared to the spectrum at maximum is the development of the 
He I 6678A and 7065A troughs. The He I reference line has 
Vco«/(He I)= 2000 km s“' and T(He I)=5 at day 19, while at day 
30 the corresponding He I parameters are Vconf(He I)= 2000 km 
s ' and T(He I)=3.4. 

Similarly to what is seen in the case of SN 1984L, the He 
I profiles retain rounded emission components that cannot be 
matched by a power-law SSp. We tried to improve the He I fits 
by switching to the e-folding assumption for the optical depth, 
introducing a continuous two-component behaviour of the He 
I optical depth: one above Vmin- 8000 km s“* with Vg =3000 
km s“' and a second component with negative {Ve - -2000 
km s '), such that t is continuous at the detachment velocity 
(i.e. at 8000 km s * for the 19day spectrum and 7000 km s * 
for the 30days one). It is important to recall here that the line 
should not be said to be “detached”. It has however a maximum 
value of T that is not at the photosphere as is ordinarily the 
case for an undetached line. The bottom panels in Figures 24 


and 25 illustrate the resulting synthetic spectra in the e-folding 
cases. The fit is somewhat improved compared to the power- 
law case. The noticeable improvements concern in particular 
the He I features and the infra-red Ca II profile. We note the 
good match to the ~6250A trough with Ha. The best fit is 
achieved using Vcon((Ha)= 7000 km s“* and T(Ha)=2.5 at day 
19 (Fig. 24), while for day 30 Ha has Vcon/(Ha)= 8000 km s“* 
and T(Ha)=1.45 (Fig. 25). No alternative identification to Ha, 
that would be logically acceptable, has been found. In addition, 
the Ha identification is strongly supported by the presence of 
the absorption notch near 4660A, well accounted for by H/? in 
our resulting synthetic spectra. 

3.2. Type Ic SNe: the representatives 
* SN 1987M: 

In Figure 26 the observed spectrum of SN Ic 1987M, near 
day 20, is compared to the generated SSp with Vpbot - 7000 
km s * and Tbb =4500 K. The OI 1113A trough is clearly deep 
and is well matched with our SSp (t =5). Identifying He I lines 
in this object is problematic and was a subject of different in¬ 
vestigations (Jeffery et al. 1991; Clocchiatti et al. 1996b). In 
fact, Jeffery et al. (1991) have presented a synthetic-spectrum 
analysis around maximum and came to the conclusion that He 
I lines may be present. In addition they also claimed the pres¬ 
ence of a weak Ha trough. In our SSp-LTE approach applied to 
the 20day spectrum, we tested two possibilities for He I identi¬ 
fication, namely He I lines at high velocity (Vco„,(He I)= 10000 
km s“'; T(He I)=0.16) and a lower velocity case (veon/(He I)= 
2000 km s '; T(He I)=0.6). The upper panel in Figure 26 cor¬ 
responds to the first case, in which the strong P-Cygni profile 
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Fig. 26. SSp fit of SN 1987M compared to the observed spec¬ 
trum at 20days. The upper panel shows the high He I velocity 
case, while the bottom panel displays the low He I velocity case 
(see discussion). Conspicuous line features are shown. 

at ~5750A is due to Na I line alone. The high velocity He I 
5876A then contributes to the trough blueward of the Na I pro¬ 
file. The undetached Sc II lines are in part responsible for the 
two absorption features indicated in the plot (t(ScII)= 0.7). The 
other He I optical lines, at 6678A and 7056A, both contribute to 
the strong Fe II features (T(FeII)=22). In the low velocity case 
(Fig. 26; lower panel) a weak He I absorption at ~6500A ap¬ 
pears, which provides a better match to the observed spectrum, 
as does “He I-i-Na I D” at ~5750A. It is not simple to decide 
which fit to adopt from analyzing only one spectrum. The two 
possibilities seem plausible and both indicate the presence of 
He I lines in SN 1987M at the age of 20 days. 

We find no need to include either Si II and/or C II or Ha. 
We note however that Si I lines may help to fit the deep trough 
around 6900A, but they would introduce then various unwanted 
features in the rest of the spectrum. 

* SN 19941: 

SN 19941 is considered the best observed “normal” Type Ic 
event. Various works have presented detailed analysis of line 
identifications in this object (Wheeler et al. 1994; Filippenko 
et al. 1995; Clocchiatti et al. 1996b; Millard et al. 1999). The 
most controversial issues concerned traces of H and He in the 
early spectra. Filippenko et al. (1995) have attributed the ob¬ 
served absorption feature near 10250A to the infrared He I 
10830A, arguing for the presence of helium in the ejecta of 
SN 19941. Millard et al. 1999, however, clarified the related is¬ 
sues by means of synthetic spectra. One of the most important 
results was the incompatibility of the simultaneous fit of the 
~10250A observed feature with He I 10830A and with the op¬ 
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Fig. 27. SSp fit of SN 19941 compared to the observed spectrum 
at maximum brightness. Conspicuous line features are shown. 
The window shows the C II and Si II identification cases (see 
discussion). 


tical He I lines. Contrary to the case of SN 1999ex (see Fig. 19), 
trying to account for the observed absorption near 10250A with 
He I 10830A introduces too strong He I features in the opti¬ 
cal region (Millard et al. 1999). Even with a NLTE synthetic 
spectra analysis, Baron et al. (1999) were not able to repro¬ 
duce the ~10250A infrared feature assuming it to be the He I 
10830A line. Instead, the feature in question may be accounted 
for by lines of C I and/or Si I (Baron et al. 1999; Millard et al. 
1999). 

Eigure 27 compares the observed spectrum around maxi¬ 
mum with an SSp that has Vphot - 12000 km s“' and Tht, =7000 
K. The match is quite good. The trough near 5700A is ac¬ 
counted for by the undetached Na ID line (t = 1). The other 
prominent features are well described by lines of Ee II(t = 10), 
Ca II(t = 300), O I(t = 0.7), [O II](t = 0.2). Lines of Mg II, 
r = 1, have been introduced to help the fit blueward the IR-Ca 
II P-Cygni profile. Lines of Ti II, t = 0.8, are also considered 
to help the fit in the 4000-4500A region. 

Concerning the feature near 6180A, it was difficult at some 
epochs in the spectra presented by Millard et al. (1999), to de¬ 
cide between detached CII6580A and undetached Si II6355A. 
Moreover, at -2days and -4days, the authors obtained a better 
fit with a combination of the two. In the spectrum displayed in 
Eigure 27 we test the two possibilities, namely undetached Si II 
(t = 1.5) and detached C II {Vcom- 8000 km s“^; t=1.7x10'^). 
A closer view of the region of interest is displayed in the win¬ 
dow. We found Si II somewhat too blue to account for the ob¬ 
served feature. There is no need at this epoch to combine the 
two lines. Therefore we prefer detached C II as the more prob¬ 
able. 
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Fig. 28. SSp fit of SN 19996aq compared to the observed spec¬ 
trum at maximum light. The right and left windows show a 
zoomed view of the 6300A and 5700A regions, respectively. 
Conspicuous line features are also reported. 

* SN1996aq: 

We analyzed two photospheric spectra of SN 1996aq, 
namely one near maximum light and one at ~24days. The 
supernova has been classified as Type Ic SN near maximum 
when first observed (Nakano et al. 1996). The observed spec¬ 
trum near maximum is compared, in Figure 28, with a synthetic 
spectrum that has Vphot - 9000 km s ' and Tbh =9000 K. The 
SSp accounts for almost all the conspicuous features, namely 
IR-Ca II, Fe II, O I and Mg II. The trough near 5680A appears 
broader than usual. To account for this absorption feature we 
use a combination of undetached Na ID (r = 1) and detached 
He I {Vcont- 3000 km s *; r = 0.35). The resulting profile fits 
nicely the broad trough (left window in Fig. 28). The absorp¬ 
tion features due to the He I lines 6678A and 7065A are too 
weak to be clearly seen as indicated in Figure 28. Concerning 
the ~6300A feature, while Si II is found to be too blue, we find 
that a combination of C II (Vcom- 2500 km s“*; t = 2.2 x 10“^) 
and Ha {Vcont- 6000 km s“'; t = 0.12) provide a satisfactory 
fit. On the one hand, the weak feature near 6330A is produced 
by the minimum of the C II6580A absorption. This possibility 
is illustrated in the right window of Fig. 28. On the other hand, 
a weak absorption in the emission feature near 4560A could be 
due to Hj6 although the Ha optical depth seems too small for 
Hy6 to be plausible. 

Alternatively, if we try to fit the 6300A feature by only C 
II, then this latter should be expelled at very high velocity, even 
more than helium. Moreover, in this way we could not produce 
the notch near 6330A. We regard the “C II-i-Ha” combination 
as the more probable. However a final confirmation of this is 
still not beyond doubt. The identification of Ha in Type Ic was a 


subject of different discussions (Jeffery et al. 1991; Filippenko 
1992; Swartz et al. 1993; Wheeler et al. 1994; Branch 2003). 
Generally the identification of Ha in early spectra of some SNe 
Ic objects (exp. SN 1987M and SN 19941) has not been ac¬ 
cepted. In SN 1987M for example the measured velocity for 
Ha was suspiciously low compared to calcium and oxygen. In 
the case of SN 1996aq we find a different situation, since the 
expansion velocity attributed to Ha is the highest (i.e. higher 
than C II and He I while the other lines are all undetached). It 
would be interesting to have more detail for this object and to 
have a larger sample of early spectra of Type Ic objects. 

Figure 29 compares the 24day spectrum to the SSp that has 
Vphot - 4600 km s“' and T;,* =5000 K. Ions that are responsi¬ 
ble for the most conspicuous supernova absorption features are 
indicated. The match is good. Note both the strong Ca II in¬ 
frared triplet and Ca II H&K are well produced in the SSp. Ni 
II helps the fit longward of the Ca II H&K profile. The double 
absorption feature near 9000A is well produced by undetached 
lines of C I(r = 1) and Mg II(t = 3). Lines of O I and [O 
II](at 7321A) are also introduced to account for features near 
7640A and 7220A, respectively. A noticeable feature of the SN 
at this phase is the clear appearance of He I lines, that question 
the classification of SN 1996aq as Type Ic object. The broad P- 
Cygni profile with the absorption minimum near 5720A cannot 
be accounted for by Na ID alone. We favour He I 5876A be¬ 
cause it implies other observed features, especially the He I 
lines at 6678A and 7065A. He I lines in the SSp are produced 
adopting Vcom- 2400 km s“* and r = 0.25. 

The slope in the 5950-6450A wavelength range does not 
contain conspicuous features although it has some weak ab¬ 
sorptions, reminiscent of what is seen for example in SNe 
1984L and 1998T at similar phases. We regard SN 1996aq as a 
transition object between Ib and Ic, rather than a "pure” Type 
Ic event. We have checked different line combinations in order 
to decide which ions to introduce to fit the 5950-6450A region, 
and especially the feature near 6355A. The different panels in 
Figure 29 illustrate these possibilities. In fact, while lines of 
Fe II and Si Il(undetached) provide a good fit to the observed 
features, we find it difficult to decide between C II, Ha and 
Ca I to account for the absorption near 6355A. The three cases 
correspond respectively to: C II with Vco«/= 5400 km s * and 
r = 7 X 10“® (upper panel); Ha with Vcont- 4900 km s“' and 
T - 0.02 (middle panel); undetached Ca I with t = 6 (lower 
panel). The Ca I case, lower panel, produces an unwanted ab¬ 
sorption trough near 6060A. For the C II case, upper panel, the 
element is assigned a velocity 3000 km s ' greater than the one 
of helium, which is hard to accept. Therefore the most prob¬ 
able situation remains the case of the high velocity hydrogen 
(middle panel). 

3.3. Type lib SNe: the representative 
* SN 1993J: 

The discovery of the "hybrid” Type Ilb SN 1993J has cre¬ 
ated a link between SN Ib-c and Type II objects, creating and 
extending our understanding of the physics of core-collapse. At 
early phases, this SN displays conspicuous hydrogen Balmer 
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Fig. 29. SSp fit of SN 1996aq compared to the observed spec¬ 
trum around day 24. The different panels show the possi¬ 
ble line identihcations for the 6300A absorption, namely C 
II (upper panel), Ha (middle panel) and Ca I (bottom panel). 
Conspicuous line features are reported. 

features similar to Type II SNe. At the nebular phase however, 
the spectrum shows many signatures of Type Ib-c SNe, namely 
strong oxygen and calcium lines (Filippenko et al. 1994; Lewis 
et al. 1994; Matheson et al. 2000). Furthermore, the analysis of 
late epoch spectra reveals that traces of hydrogen (i.e. Ha in 
emission were still present (Patat et al. 1995). Another well ob¬ 
served SN belonging to Type Ilb class is SN 1996cb. By com¬ 
parison with SN 1993J, SN 1996cb showed Balmer lines with 
stronger P-Cygni prohles. In addition, the photosphere of SN 
1996cb receded faster than its counterpart in SN 1993J and the 
He I features (in absorption) appeared earlier, around day 24 in 
SN 1993J and near maximum light in SN 1996cb (Qiu et al. 
1999; Deng et al. 2001). 

We have analyzed spectra at three different epochs for SN 
1993J, namely at 16, 24 and 59 days. Figure 30 displays our 
best ht spectra compared to the observe ones. The plotted syn¬ 
thetic spectra have VpUot - 9000 km s“' and T** =7800 K 
(16days; upper panel), Vphot - 8000 km s“^ and Tbb =7000 K 
(24days; middle panel) and Vphot - 6000 km s“' and Tbb =5000 
K (59days; lower panel). The ht to the 16day spectrum is 
good, except the Ha prohle. The corresponding strong and 
broad P-Cygni feature cannot be produced completely by the 
“SYNOW” code. This is because in our SSp treatment we are 
adopting a resonance scattering source function. The hydrogen 
reference line has Vcom- 1000 km s ', while an optical depth 
of T = 21 is used. Prohles of H/?, Hy and Hh are clearly dis¬ 
cernible. The P-Cygni prohle around 5700A is assigned to He 
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Fig. 30. SSp ht of SN 1993J compared to the observed spectra 
at 16days (upper panel), at 24days (middle panel) and at 59days 
(lower panel). Conspicuous line features are reported. 

15876A, with a small contribution from Na ID. The other opti¬ 
cal He I lines are not visible at this phase. The helium is found 
to be undetached and has a moderate optical depth of 0.6. The 
shallow IR-Ca II triplet is htted adopting r = 20. The observed 
Ca II H&K is nicely reproduced in the SSp. 

As the supernova evolves, the IR-Ca II prohle becomes 
clearer (Fig. 30; middle panel). At this phase an optical depth of 
T(CaII) - 300 is used. Fe II lines are also clearer compared to 
the 16day spectrum (T(Fell) = 3 at 16days and T(Fell) =12 
at 24days). The notch appearing on the emission component 
of Ha is attributed to He I 6678A. Helium is still undetached 
with an optical depth of 1.5. SN 1996cb at a similar phase, near 
25days, has more evident He I lines compared to SN 1993J, es¬ 
pecially the lines at 6678 and 7056 (Fig 2; Deng et al. 2001). In 
fact, its 25day spectrum resembles the 59day spectrum of SN 
1993J rather than the 24days one. The He I lines at day 59 are 
prominent (Fig. 30; bottom panel). The He I reference line is 
undetached and has T(Hel) = 30. The corresponding SSp ht is 
quite good for the 5876A line, while it is difficult to obtain a 
good match with the other optical He I lines. 

Hydrogen Balmer lines become weak at this phase. Ba II 
lines are introduced to improve the ht just blueward of He I 
5876A (T{BaII) = 5). At this phase, forbidden emission lines 
develop, indicating the transition to the nebular phase (exp. 
[O I] 5577A , [O I] 6300,64A and [Ca II] 7291,7324A; bot¬ 
tom panel). The O I 7773A absorption prohle is reproduced at 
the three different epochs by t = 0.3(at 16days), t = 0.4(at 
24days) and t = 2.5(59days). Ni II lines are also introduced in 
the 24days and 59days synthetic spectra, blended with Hd, to 
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improve the match with the observed weak absorption on the 
Ca II H&K emission component. 

3.4. Type II SNe: the representative 

* SN 1999em; 

As a representative of the Type II class, we have analyzed 
SN 1999em. The well sampled observations, spectra and pho¬ 
tometry, have been presented and studied by Elmhamdi et al. 
(2003). The event shows characteristics typical of a Type II-P 
object, namely clear and broad Balmer P-Cygni profiles and a 
plateau phase of almost constant luminosity in the optical. 

Figure 31 displays two spectra during the photospheric 
phase. For comparison with other events of our CCSNe sam¬ 
ple, the phases are normalized to maximum light in Type Ib-c 
SNe, adopting 16 days as the rise time to reach maximum. The 
observed spectra are compared to the best fit synthetic ones that 
have Vphot - 10000 km s“' and Tht, =10000 K (-6days; upper 
panel) and Vpi,gt - 4600 km s“^ and T^^, =6400 K (25days; 
lower panel). Synthetic line features are designated according 
to the ion whose line gives rise to the feature (Fig. 31). At the 
earliest phase, i.e. -6days, only three elements are introduced 
in the SSp model. Indeed undetached lines of Balmer hydro¬ 
gen, He I and a weak contribution from Ca II are almost suffi¬ 
cient to reproduce the most conspicuous features superimposed 
on the “hot” continuum. The fit with Na I D, at this phase, is 
poor compared to He I 5876A. Moreover He I contribution, in 
Type II SNe, is found to be important shortly after the explo¬ 
sion. In fact in SN 1987A, the He I 5876A feature was clearly 
present during the first few days, then owing to the decreasing 
envelope-temperature it rapidly faded and disappeared com¬ 
pletely around 1 week after the explosion, when Na I D starts 
to emerge (Hanuschik & Dachs, 1988). Baron et al. (2000), in 
analyzing very early spectra of SN 1999em, have found evi¬ 
dence for helium enhanced by at least a factor of 2 over the 
solar value. The interpretation of the He I strength in terms 
of helium-overabundance is however premature, since freeze- 
out effects were not considered and could lead to the enhanced 
helium-excitation compared to the steady-state model (Utrobin 
& Chugai 2002). 

Balmer hydrogen features, with r = 15, are clearly evident. 
As in the case of SN 1993J, we cannot reproduce the full ob¬ 
served Ha profile for the same reason. At day 25, Ha P-Cygni 
profile becomes narrower because of the decrease in the ex¬ 
pansion velocity (Fig. 31; bottom panel). Ha in the SSp has 
now Veo„f= 1000 km s ' and t = 21. The envelope tempera¬ 
ture decreases and many lines emerge at this phase. Apart from 
the hydrogen lines (slightly detached), all the lines are unde¬ 
tached. The match with almost all observed features is good. 
The Ca II H&K and IR-Ca II are nicely reproduced. A few dis¬ 
tinct absorption lines redward of the near infrared Ca II triplet 
are identified as C I, t = 0.6, through our spectral synthesis. 
Sc II lines, r = 3, account rather well for features near 5450, 
5600 and 6170A. The blue part of the spectrum, A < 5000A, 
is subject to severe line blending, but it seems well matched 
with contributions from lines of Fe II(t = 10), Ni II(t = 5), Ti 
II(t = 8), Ba II(t = 2), Mg II(t = 2), Ca II and hydrogen. The 
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Fig. 31. SSp fit of SN 1999em compared to the observed 
spectra at day -6 (upper panel) and at day 25 (lower panel). 
Conspicuous line features are shown. For comparison, the re¬ 
ported phases are normalized to maximum in Type Ib-c SNe, 
using “16 days” as the rise time. 

feature near 7670A is attributed to O I 7773 A(t = 0.3), while 
there may be a contribution from [O II] around 7230A. 

4. Spectroscopic mass estimates 

In the following we adopt two methods in order to obtain 
rough estimates of some ejecta masses (i.e. the total ejecta mass 
and/or a given element mass). The methods use the results from 
our spectral best-fits. 

4.1. Method 1 

We compare the optical depths derived from our best-fits to the 
optical depth plots for different compositions presented in the 
work by Hatano et al.(1999), making sure of course that the 
same reference lines are used. In that paper the authors pre¬ 
sented a systematic survey of ions that could be responsible for 
supernova features in six different compositions. The FTE op¬ 
tical depth of each reference line, of a given ion, is computed 
and then plotted against temperature. We focus on the ratio of 
hydrogen to helium optical depths, i.e t(H I)/T(He I), attempt¬ 
ing thereafter to translate this ratio into a relative abundance of 
the two elements. This is important as it can give us a rough 
idea about the amount of hydrogen present in these objects. 
The optical depth ratio is assumed to trace the abundance ra¬ 
tio of the two elements. This is essentially due to the fact that 
the T-calculations of Hatano et al. (1999) adopted a fixed elec¬ 
tron density. The r(H I)/T(He I) ratio, indeed, depends on elec¬ 
tron density, temperature and hydrogen to helium abundance 
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ratio. Therefore for a given temperature, the t(H I)/T(He I) ra¬ 
tio depends in principle only on the abundance ratio. We must 
remember however that the “r-plots” are for LTE, with no al¬ 
lowance for non-thermal excitation is made. For helium, for 
example, non-thermal excitation effects are significant in de¬ 
veloping He I features (Lucy 1991). More detailed estimates 
might take this effect into account, if one has a clearer under¬ 
standing of outward mixing of radioactive material. 

For SN 19901 around maximum and according to our best 
SSp fit, the optical depth ratio, evaluated at the same velocity^, 
is measured to be r(H l)/T(He I)~0.6. At similar high tempera¬ 
tures (i.e. ~ 14000 K), one obtains a value of about 40 for that 
ratio according to the optical depth plots of different ion ref¬ 
erence lines for a hydrogen-rich composition in Hatano et al. 
(1999) (see their Fig. 2). In addition, in the cited work the t(H 
I)/T(He I) ratio remains almost constant at high temperatures, 
then it starts increasing once the temperature decreases. For ex¬ 
ample around T~ 10000 K, the ratio is about 80. Near 9000 K, 
however the ratio reaches a value of about 1000. These val¬ 
ues (i.e. 40, 80 and 1000 at their corresponding temperatures) 
correspond to a factor of 10.23 in abundances between hydro¬ 
gen and helium (Table 1; Hatano et al. 1999). Assuming pro¬ 
portionality, we obtain then a rough estimate of the hydrogen 
to helium relative-abundance for the case of SN 19901 around 
maximum. In fact a relative abundance of ~0.15 is found. For 
SN 2000H in which Her is clearly deep (t~ 5 at maximum; Fig. 
23), an abundance ratio (i.e. H/He) is estimated to be ~0.23. 
SN Ilb 1993J, around day 16, has a ratio ~0.8. We estimate 
an even higher ratio for SN 1999em, namely ~3.2. These esti¬ 
mates, although rough, are in accord with our impression that 
the hydrogen mass increases as we go from Type Ib to lib to 
Type II SNe. Figure 32 displays the logarithmic hydrogen to 
helium abundance ratios against temperature for the events of 
our sample. Note here that for low temperatures, T < 9000, 
the reported values in Figure 32 can be taken as upper lim¬ 
its. This is because we used a value of 1000 for the ratio t(H 
I)/T(He I). At low temperatures, however, the ratio becomes 
significantly larger than this adopted value (Fig. 2a; Hatano et 
al. 1999). An additional source of uncertainty may be related 
to the continuum temperature estimates owing to the total red¬ 
dening effects. Furthermore, according to our fit experience on 
the sample spectra, an uncertainty of about 10 % is assigned to 
the derived values of the optical depths. In the case of SN IIP 
1999em and SN lib 1993J, with their complete Ha P-Cygni 
profiles, greater optical depth uncertainties should be expected, 
although in these events we got good fits to the H/1 profiles. 
It is worth noting however that the enormous NLTE effects in 
the He I lines lead to very inaccurate determination of helium 
abundances. 

Making use of the derived photospheric velocities “Vpho", 
from our synthetic spectra modeling, it is possible to recover 


^ The optical depths reported in Table 2 are evaluated at the detach¬ 
ment velocity of the corresponding reference lines (e.g. for SN 19901 
at maximum t(H I) at 16000 km s^'and T(He I) at 14000 km s“‘). For 
a “v^*” power-law, the optical depth of He I evaluated at the same ve¬ 
locity of HI is then: 
r(HeI)= 2.9x(14/16)* 1.0 
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Fig. 32. The hydrogen to helium abundances ratio against tem¬ 
perature. See Sect. 4 for more details 


“spectroscopic” estimates of the kinetic energy and mass 
above the photosphere. In fact, Millard et al. (1999) have shown 
that for spherical symmetry and an“r density distribution, 
the mass (in Mq) and energy (in 10 ^^ ergs) above the electron¬ 
scattering optical depth “ Tg/’ can be expressed as: 


M = (1.2 X 10 '^) V 4 Pe Tes fmin) 

( 1 ) 

E^(\.2x\0-^)v\tlpeTesfE{ny, 

( 2 ) 


where fM(n) - ^ and /^(n) = td is time after explosion 
in days, V 4 is Vpi,ot in units of 10 ^ km s“' and is the mean 
molecular weight per free electron. 

In applying the above equations to our supernova sam¬ 
ple, we adopt n = 8 and Tes - 2/3. For the mean molecular 
weight per free electron we have different values for each case, 
namely: SNe Ib: Pe - ^ for half ionized helium or doubly ion¬ 
ized oxygen; SNe Ic: pe - 14 for a mixture of carbon and oxy¬ 
gen, both singly ionized; SNe II: p^ - \ Xo 2 for fully ionized 
or half ionized hydrogen. For the case of “hybrid” lib SNe we 
assume a mixture of hydrogen and helium, both half ionized 
and hence pe - 5. 

Figure 33 displays results for five events. Filled symbols 
refer to the mass (in units of M©; upper panel), while the open 
ones indicate the kinetic energy (in units of 10 ^* ergs; lower 
panel). For each supernova, the estimated masses for different 
epochs are connected by a dotted line. The short-dashed line 
connects the derived energies in order to clarify trends. For a 
given phase, the reported amounts indicate the mass moving 
above the photospheric velocity, carrying a corresponding ki¬ 
netic energy. The events displayed in Fig. 33 have been selected 
on the basis of their type and velocity. SN 19901, being an ob¬ 
ject with high velocities, appears to have high kinetic energy, 
while SN 199ID with its lower velocity structure lies at the bot¬ 
tom of the plot. For a significant comparison, we choose data at 
similar phases, around day 38. At that epoch, the mass moving 
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Days since explosion 


Fig. 33. The evolution in time of the “spectroscopic” estimates 
of Mass (upper panel) and Energy (lower panel) above the pho¬ 
tosphere. 

faster than “Vphot' and the corresponding energy are estimated 
to be: (1.1 Mo; 1.6 foe) for SN 19901; (0.45 M^, 0.28 foe) for 
SN 2000H; (0.3 Mq; 0.1 foe) for SN 1991D; (0.6 Mq; 0.8 foe) 
for SN 1993J and (0.95 foe) for SN 19941. SN lb 

1984L has behaviour similar to SN 2000H. These derived 
quantities suffer, of course, from some uncertainties, namely 
those incurred with the derived photospheric velocities'^ from 
our best fits and the adopted value for the power-index “n”. 
However, they seem to give reasonable and meaningful val¬ 
ues that are not in disagreement with other methods (e.g. light 
curve modeling ; NLTE treatment of early and late spectra), 
especially for the well studied objects such as SNe 19901, 
19941 and 1993J. The oxygen-mass estimates could be a further 
check of the trend in the derived photospheric outflow mass. 
Elmhamdi et al. (2004), on the basis of [O 1] 6300, 6364A line 
analysis at nebular phases, have estimated a lower limit on the 
oxygen mass to fall in the range 0.7-1.35 Mq. The available 
estimated amounts for the other objects are as follows: ~0.3 
Mo for SN 1984L (Eilippenko et al. 1990); -0.5 Mq for SN 
1993J (Houck & Eransson 1996) and -0.4 Mq for SN 19941 
(Woosley et al. 1995). These estimates follow the trend seen in 
the masses above the photosphere. 

We emphasize here that this simple approach, using equa¬ 
tions 1 and 2, does not provide “reasonable” estimates in the 
case of Type 11 SNe. In fact when this method is used for both 
SNe 1987A and 1999em, the estimated “M” and above 
the photosphere seem to be very small, in contradiction with 

^ A mean fitting uncertainty of -500 km s“' is assigned to the de¬ 
rived velocities. The uncertainty in the velocity would affect more the 
uncertainties in the derived energies than those of the derived masses. 


what we may expect. Two possible effects might invalidate the 
method for Type 11 objects. Eirst, because the method cannot 
be used whenever the power index “n” value falls below 5 
(i.e. to keep /«(«) and /^(n) positive and convergent). Eor SN 
1987A, indeed, Jeffery & Branch (1990) had found “n” falling 
below 5 already within a week after the explosion (see their 
Eigures 15 and 19). Second, in Type 11 SNe and after the earliest 
times, the electron-scattering optical depth for the photosphere 
comes from a thin layer of ionized hydrogen at the recombina¬ 
tion front. Earther out, there can be a lot of neutral hydrogen 
(and helium) that is not accounted for in the “Millard et al.” 
equations adopted in our simple analysis. When, for instance, 
“/le” is set to one, the equations are only giving the mass and 
energy of the matter near the photosphere. Consequently one 
may question why they should be taken seriously for the other 
CCSNe classes. Some carbon and oxygen could be neutral too. 
One simple explanation could be that because of the large hy¬ 
drogen envelope, these effects seem to be more “severe” in 
Type 11 events rather than the rest of SNe. In summary, the sim¬ 
ple method does not apply to SNe 11 because it does not take 
into account the recombined hydrogen. 

By adopting the estimate of the outflow mass above “Vphot” 
as an upper limit of the helium mass in Type lb and 11b events, 
and using the derived constraints on the hydrogen to helium 
abundances ratio, one can thus recover an upper and rough 
estimate of the mass of expelled hydrogen. Around maxi¬ 
mum light, the hydrogen amount in SN 19901 is found to be 
-0.16 Mq. At similar phase, an amount of about 0.1 M© is 
computed for SN 2000H. SN 1993J, however, seems to eject 
a quantity as high as -0.7 Mq. SN 1993J indeed serves as a 
control case for our adopted methodology. Different studies ar¬ 
gued that a low-mass hydrogen envelope on a helium core is 
the most likely scenario for the progenitor of SN 1993J. On 
the one hand, light curve modeling demonstrates a compati¬ 
bility fit with a 4 Mq helium core and hydrogen envelope of 
-0.2 Mq (Woosley et al. 1994). A similar analysis invoked 
a helium core -4-6 Mq, with a residual hydrogen mass less 
than 0.9 Mq (Nomoto et al. 1993; Shigeyama et al. 1994). On 
the other hand, best fit models of spectra lead to similar pro¬ 
genitor properties, namely 3.2 Mq for the helium core with a 
0.2-0.4 Mq hydrogen envelope (Houck & Eransson 1996). Our 
estimate is therefore in good agreement with the previous cited 
studies. Eor the rest of Type lb objects, an upper limit of the 
hydrogen mass of the order 0.1 Mq is estimated. 

4.2. Method 2 

Alternatively, an approximate way to estimate the hydrogen 
mass invokes the amount needed to fill a uniform density 
sphere of radius “vxf” at an epoch since explosion. Adopting 
a value of “10^®” for the hydrogen Balmer fraction (within the 
LTE context), the required ion mass can be given by: 

M(Mo) - (2.38 X 10-^) tj T{Ha)- (3) 

This comes from the equation for the Sobolev optical depth for 
an expanding envelope (see e.g. Castor 1970; Jeffery & Branch 
1990). where td in days and va, is in lO'^ km s“*. This equation is 
then used to determine the amount filling a spherical shell of an 
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inner radius “v„„„ x t/’ and an outer edge of radius ‘"v^ax x td”- 
Constraints on the shell widths, - v„,ax”, are approximated 
using Gaussian fits to the observed Ha absorption troughs by 
means of the corresponding FWHM. It is worth noting here that 
constraints on “v,„ax” from this method, i.e. Gaussian fit, can be 
considered as a lower limit to the “real” maximum velocity. 

Although non-thermal excitation and “NLTE” effects may 
be also important for hydrogen, the method seems to give rea¬ 
sonable estimates. For SN 19901 around maximum light a value 
of 0.02 Mq is computed, while at similar phase SNe 1983N and 
2000H have, respectively, 0.008 and 0.08 Mq. For type Ilb SN 
1993J, at day 16 since maximum, the method gives an approx¬ 
imate amount as high as 1.7 Mq. This amount is of course very 
sensitive to the assumed optical depth. In fact SN 1993J has an 
optical depth similar to Type IIP 1999em (i.e. T{Ha) - 21 at 
~25days). However the fit is more “convincing” in SN 1999em 
than in SN 1993J, especially in the Ha absorption troughs (see 
Fig. 30-middle panel and Fig. 31-lower panel). This might be 
indicative of an overestimated optical depth in SN 1993J and/or 
a larger departure from “FTF” compared to SN 1999em. 

A rise time of 15 - 20 days has been adopted in this class of 
object. Moreover, if we adopt a representative Ha optical depth 
of 0.5 at day 20 with an Ha velocity restriction similar to SN 
19901, the estimated hydrogen mass is of the order 0.015 Mq. 
Events with higher velocity widths and/or deeper Ha troughs 
would eject larger amounts. 

5. Discussion and conciusion 

One of our main goals was to identify traces of hydrogen in 
Type Ib-c, and to identify any systematic similarities and differ¬ 
ences correlating with other physical properties. Consequently, 
guided by modeling early spectra of SN Ib 19901, we have ex¬ 
plored different combinations of ions and shaping the 6000 - 
6500A wavelength range. Special attention has been devoted to 
the feature seen near 6300A. Even though the number of events 
with spectra well sampled in wavelength, is limited, it appears 
that hydrogen in varying amounts is identifiable in most Type 
Ib objects especially at very early times. Only in two cases, 
namely SNe 1991D and 1991F, does the Ne I 6402A line re¬ 
main as an alternative possibility with large departure from 
FTF. Even in SN 1996aq, that we classified as a transition Type 
Ib/c event, and SN 1999ex the presence of the Ha trough seems 
highly prefetTed. 

At later phases, more than ~ 3 weeks, SNe 1984F, 1988F, 
1991ar, 1998dt, 1998T and 1999dn behave quite similarly, 
namely still showing evident He I lines and a “flat” 6000 - 
6500A region of the spectrum, with a weak Ha absorption fea¬ 
ture. Even SN 199ID with its particularly low photospheric ve¬ 
locity and narrow lines belongs to this class. SN 1991F appears 
to be similar to SN 199ID in having shallow He I optical fea¬ 
tures. SNe 1999di and 2000H are the unique events among the 
sample that display obvious He I lines, together with a pro¬ 
nounced and deep Ha absorption line. It is interesting to note 
here the peculiar behaviour of SN 1997dc. The event has He I 
troughs as deep as in SNe 1999di and 2000H, however it has 
a flat 6000 - 6500A region. SN 19901 at this phase shows dis¬ 
tinct He I lines and a moderate Ha trough. Figure 34 compares 



Wavelength (A) 


Fig. 34. Comparison of SN 19901 with SNe 1984F and 2000H, 
at two different phases of evolution. The vertical dashed lines 
indicate the locations of Ha and H/l troughs in SN 2000H. He 
I optical series are marked by vertical ticks. 


SN 19901 with SNe 1984F and 2000H at two different phases 
of evolution. The vertical dashed lines indicate the locations of 
Ha and H/? troughs in SN 2000H. On the one hand, while at 
early phases (upper panel) the trough assigned to Ha is clear 
in all the three SNe, it disappears later on in SN 1984F (lower 
panel). On the other hand, appears clearly in SN 2000H 
and absent in SNe 1984F and 19901. Moreover, the comparison 
emphasizes the high velocity behaviour of SN 19901 as is evi¬ 
dent from the blueshifted features in SN 19901 compared to the 
others. He I optical lines are also shown. For the three events 
the He I troughs at 5876A, and 7065A are clearly identified 
and get narrower with time. At early phase the He I 6678A is 
recognized in SNe 19901 and 1984F while it is just hinted in 
SN 2000H. Eater on. Fig.34-bottom panel, the He I 6678A line 
grows in strength in SNe 1984F and 2000H while it fades in 
SN 19901. The comparison of these three events, in terms of 
changes in line visibility and velocity, demonstrates the com¬ 
plexity involved in demonstrating any thing like a continuous 
sequence. 

Hydrogen manifests its presence in a different way in the 
other CCSNe types. A part from SN 1996aq that we re-classify 
as a transient Ib/c event rather than a pure Type Ic SN, there 
is no evidence for the presence of hydrogen in Type Ic ob¬ 
jects analyzed here (i.e. SNe 1987M and 19941). The “hybrid” 
SN lib 1993J displays typical Type II features at early phases, 
such as strong and broad Ha P-Cygni emission component, 
which is absent in Type Ib objects. This can be explained within 
the context of the “detachment” concept. In fact, Ha P-Cygni 
profile would lose its obvious emission component when it is 
highly detached. At early epochs, both SN lib 1993J and SN II 
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Fig. 35. The photospheric velocity evolution of Ib SNe sam¬ 
ple, compared to SNe 1987A, IIP 1999em and Ic 19941 (upper 
panel). Middle panel; the Ha contrast velocity ratio evolution 
(i.e. v^“"'’(Ha)/Vp/,o,). Lower panel: the Ha contrast velocity 
evolution (i.e. Vm;„(Ha)-Vy,/,o,). Results for SN IIP 1999em and 
SN Ilb 1993J are also shown. 


1999em show clear evidence of the other HI Balmer lines since 
the corresponding optical depths are too large to allow them to 
be distinctly visible, contrary to normal Type Ib objects, with 
the exception of some cases with deep and conspicuous Ha 
troughs that present signature of H/3 too (e.g. SNe 1999di and 
2000H). Two factors make Hy6 barely discernible in Type Ib; 
the optical depth found to fit Ha is small and the contrast ve¬ 
locity of Ha is high. The opposite is seen in Type II and Ilb 
events. As time goes on, the Ha emission peak in SN 1993J 
changes to a double-peaked structure which we recognized as 
the emergence of the He I 6678A, with more clearer He I op¬ 
tical lines. These laters disappear at the nebular phase and the 
spectrum is dominated by typical Ib features such as [O I], [Ca 
II] and Ca II in emission. The situation in SN 1999em, and in 
Type II in general, is different. The He I lines are found to dis¬ 
appear very early, about one weak after explosion, and even at 
later nebular epochs the Ha is still the most prominent feature 
in the spectrum. 

In Figure 35, upper panel, we report the resulting photo- 
spheric velocities from our best fits for the CCSNe sample. 
Data for SN 1987A, corresponding to Fe II 5018A absorp¬ 
tion, are also shown for comparison (dashed line; Phillips et al. 
1988). Additional points for SN 19941 are taken from Millard 
et al. 1999. The plot indicates the low velocity behaviour of 
Type II SNe, both 1987A and 1999em, at early as well as in- 
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Fig. 36. The Ha velocity contrast evolution, v™],“(Ha), versus 
Ha optical depth. Results for SNe IIP 1999em and Ilb 1993J 
are displayed for comparison. Note that results are displayed 
independently of the phase. 


termediate epochs. SN lib 1993J follows somewhat similar be¬ 
haviour as SN Ic 19941 in having higher velocities. As far as 
Type Ib SNe are concerned, they appear to display a different 
velocity evolution. The scatter seems to increase at intermedi¬ 
ate phases (around 20-30 days). This fact can be simply due to 
the paucity of available observations outside that range. Around 
day 20, for example, a scatter as high as 5000 km s“^ is mea¬ 
sured. SNe 19901 and 1998dt belong to a class with the higher 
“vphor”, while objects such as SNe 199ID and 1996aq have the 
lowest estimated velocities, approaching even Type II objects. 
The remainder of the Ib events follow a similar trend, namely 
the one described by Branch et al. (2002). SN 1999ex is found 
to belong to this class. Surely more data are needed to obtain 
meaningful statitical conclusions. 

The middle and bottom panels in Fig. 35 display the evo¬ 
lution the of Ha contrast velocities, and “vcont” respec¬ 

tively, for the sample. As discussed before, when discussing 
differences in the Ha profile in CCSNe, we would expect an 
increasing value of “veo„,(Ha)” going from Type II to Ilb to Ib 
SNe. This trend is in fact illustrated in the plot (bottom panel). 
A similar trend is seen in the “v™J]“(Ha)” evolution as indicated 
by the middle panel. While in SN Ilb 1993J and SN IIP 1999em 
the line is found to be either undetached or slightly detached, it 
is highly detached in Type Ib events. Moreover, the “vconf(Ha)” 
is found to increase within the first 15 days, reaching values 
as high as 8000 km s ', and then follows an almost constant 
evolution. According to Table 2 and up to ~60 days. Type Ib 
SNe have hydrogen down to 11000-12000 km s \ while in SN 
1993J hydrogen is down to 8000 km s"'. SNe II appear to have 
hydrogen down to even lower velocities (~5000 km s ' in SN 
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Fig. 37. The optical depth of the O I 7773A line according to 
the best fit spectra. Data for different SN types are reported for 
comparison. 


1999em). In addition, hydrogen in Type Ib SNe is found to have 
very small optical depths independently of the contrast velocity 
(also independently of the phase). This is shown by Figure 36. 
In fact. Type Ib objects populate the region of extremely low 
optical depths, although the “v™j,“(Ha)” spans a large range, 
contrary to SNe 1993J and 1999em. 

Table 2 also displays properties of the helium line. In Type 
Ib SNe, He I lines are found to be not always detached as is 
the case for hydrogen. The He I 5876A line is found to be 
clearly distinguishable in Type Ib objects. A contribution from 
Na ID is argued for in some cases, such as the broad feature in 
SN 1996aq and cases with an observed bump shortward of the 
minimum absorption as seen in SNe 1990B and 1997dc. It is 
important to recall here that in some cases we obtained a better 
fit to the He I optical lines adopting an e-folding SSp, which 
allows a two-component treatment of the line and a gradually 
decreasing optical depth below the detachment velocity, rather 
than a discontinuity in it. In the case of SN 1999ex, we extend 
our He I line fits beyond 1 micron, confirming the presence of 
various IR lines, namely at 1.083, 1.284, 1.700(7), and 2.085 
microns (Fig. 20), in accord with optical He I parameters. In 
our analysis, we also pointed to some interesting line identifi¬ 
cations. For example the contribution of Sc II lines in nearly 
all the CCSNe spectra, especially for reproducing the double 
absorption features blueward of the He I 5876A line. We note 
here that except for He I and H I, almost all the lines in the 
CCSNe synthetic spectra are undetached. 

Another important study is the behaviour of the O I 
7773A line. Figure 38 reports the resulting optical depths cor¬ 
responding to our best synthetic spectra fits. It is not easy to 
draw clear-cut conclusions from this figure since one needs. 


for instance, to populate the figure with more Type Ic objects. 
However, at intermediate phases, it seems that Type Ib objects 
tend to concentrate in the low optical depth region, while SN 
Ic 1987M is found to display the deepest profile. SNe 1993J 
and 1999em, at similar phases, are the objects with the low¬ 
est O I 7773A optical depth. At somewhat later epochs, tran¬ 
sient Type Ib/c objects display deep O I 7773A troughs. The 
stronger and deeper permitted oxygen lines at early phases of 
SNe Ic and Ib/c spectra might imply that they are less diluted 
by the presence of a helium envelope. Indeed one might expect 
oxygen lines to be more prominent for a “naked” C/O progen¬ 
itor core. Despite the paucity of well sampled CCSNe obser¬ 
vations, two observational aspects tend to reinforce this belief; 
First, the forbidden lines, especially [OI]6300, 6364A, seem 
to appear earlier following a SNe sequence “Ic-Ib-IIb-II”. In 
fact the oxygen line emerges at an age of 1-2 months in Type 
Ic SN 1987M (Filippenko 1997). SN Ic 19941 displayed evi¬ 
dence for the line at an age of 50 days, although some hints 
may even be seen in the ~36day spectrum (Clocchiatti et al. 
1996b). While in SN Ib 19901 it was hinted at the 70day spec¬ 
trum (Elmhamdi et al. 2004). In other Type Ib SNe it appears 
earlier than in SN 19901. In SN Ilb 1993J, a transition object, 
the line was visible in the 62day spectrum (Barbon et al. 1995). 
SN 1996cb, another well observed lib event, showed evidence 
of the [OI]6300, 6364A line around day 80 (Qiu et al. 1999). 
In SNe II, however, the line appears later: around day 150 in 
SN 1987A (Catchpole et al. 1988) and after day 138 in SN 
1992H (Clocchiatti et al. 1996a). In SN II 1999em it is sug¬ 
gested at a somewhat earlier phase compared to SNe 1987A 
and 1992H, namely at day 114. Whether all this can be un¬ 
derstood in terms of the lower progenitor mass of SN 1999em 
and its presumed lower oxygen mass remains to be investi¬ 
gated (Elmhamdi et al. 2003). Second, it seems that the neb¬ 
ular emission line decreases in breadth following the SNe se¬ 
quence above. Of course much work is still needed in this re¬ 
spect based on larger CCSNe samples at both photospheric and 
nebular phases, relating them to the photometry of the CCSNe 
variety (Elmhamdi & Danziger, in preparation). Eurthermore, 
one important point is to check for any possible correlations 
between the oxygen minimum velocity (at early epochs), their 
line widths, and to relate this to a progenitor mass indicator [Ca 
II]/[0 I] ratio (at nebular phases; Eransson & Chevalier 1989), 
the progenitor properties (i.e. masses, energies). Other factors 
such as mixing, variation in envelope densities and metallicity 
may play an additional complicating role. 

Einally we presented two methods to determine the ejecta 
and hydrogen mass in CCSNe and especially in Type Ib events. 
Although the methods are very approximate^, our results do not 
conflict with more detailed estimates, those based on hydrody- 
namical and NLTE models. A thin layer of hydrogen, ejected 
at high velocities down to 11000-12000 km s“\ appears to be 
present in almost all the Type Ib events studied here. These re¬ 
sults suggest possible directions for further more sophisticated 
work. The necessity to introduce lines of Ne I, Sc II, Ba II, Ca 
I, Ee II in some cases but not in others must inevitably raise 


^ “Method 2” is a direct method, while “Method 1” is undirect and 
gives upper limits on the hydrogen mass. 
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concerns about the identifications and other reasons for these 
variations. 
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Table 2. Main data and fitting parameters of the CCSNe sample spectra 


Supernova 

Host Galaxy 

'^recession 

Phases 

phot 

Tbb' 

t(H) 

T(Hel) 



Nbr(ions) 

(Type) 



(days) 


(K) 



(kms^^) 

(kms^^) 


SN19901 (Ib) 

NGC 4650A 

2902 

max 

12000 

14000 

0.6 

2.9 


■EHiran 

7 




12 

10000 

5500 

0.6 

1.85 

I 


7 




21 

9500 

5400 

0.26 

0.94 

■ mm 

i^^i 

7 

SN1983N (Ib) 

NGC 5236 

513 

max 

11000 

8000 

nK|H 


15000 

^phot 

8 




10 

7000 

5000 



12000 

9000 

11 

SN1984L (Ib) 

NGC 991 

1534 

8 


HI 



■Hiiiin 

10000 

5 




32 


■1 





5 

SN1987M- (Ic) 

NGC 2715 

1339 

20 

7000 

4500 

— 

0.6 

— 

9000 

7 

SN1988L (Ib/c) 

NGC 5480 

1856 

20 

8000 


0.7 

1 

16000 

11000 

9 

SN1990B (Ib/c) 

NGC 4568 

2320 

8 



0.5 

0.55 

13000 

10000 

7 




30 

■I 

BckiI 

— 

0.8 

— 

10000 

9 

SN1991ar (Ib) 

IC 49 

4562 

28 

8000 

EEHin 


10 

16000 

^phot 

10 




48 

7000 



4.8 

16000 

8000 

10 




60 

6000 


0.28 

4.4 

14000 

7000 

12 

SN1991D (Ib) 

PGC 84044 

12500 

21 

4600 

7000 

0.46 

1.8 

12000 

6000 

11 

SN1991L (Ib/c) 

M 7-34-134 

9050 

31 

5000 

6000 

0.34 

2 

13000 

^phot 

12 

SN1993J (Ilb) 

M81 

35 

16 


■gun 

21 


MHiiiin 

^phot 

10 




24 

8000 


25 

1.5 


^phot 

10 




59 



2 

30 


^phot 

11 

SN1994I (Ic) 

NGC 5194 

461 

max 

12000 

7000 

— 

— 

— 

— 

9 

SN1996aq (Ib/c) 

NGC 5584 

1602 

2 

Hi 


0.12 


15000 

12000 

8 




24 

IBI 



Hi 

— 

7000 

11 

SN1997dc (Ib) 

NGC 7678 

3480 

28 

5000 

4200 

0.2 

10 

13000 

7000 

10 

SN1998dt (Ib) 

NGC 945 

4580 

8 

9000 

5600 

0.3 

4 

17000 

11000 

6 




33 

9000 

5000 

0.2 

10 

17000 

^phot 

9 

SN1998T (Ib) 

NGC 3690 

3080 


HI 

mjiji 


MM 

10500 

9000 

9 






Bnal 


■iIBi 

— 

8000 

9 

SN1999di (Ib) 

NGC 776 

4920 

21 

■buiim 

4800 

3 

12.6 

12000 

8500 

8 




45 


5200 

1.6 

10 

12000 

7500 

9 




52 


5000 

1 

10 

12000 

7000 

9 

SN1999dn (Ib) 

NGC 7714 

4920 

-10 

14000 

6600 

1.5 

2 

18000 

^phot 

5 




max 

10000 

6000 

1 

5 

14000 

11000 

6 




10 

7000 

5800 

0.24 

1.9 

15000 

9000 

10 




38 

6000 

5400 

0.34 

14.5 

14000 

7000 

10 

SN1999ex (Ib) 

IC5179 

3498 

4 

10000 

5800 

1.45 

2.35 

■grain 

11000 

7 




13 

7000 

5600 

1.6 

3.4 


8000 

12 

SN1999em (IIP) 

NGC 1637 

717 

-6 

10000 

10000 

15 

0.6 

^phot 

^phot 

3 




25 

4600 

6400 

21 

— 

5600 

— 

13 

SN2000H (Ib) 

IC 454 

3894 

max 

11000 

9000 

5.25 

2 

13000 

^phot 

9 




5 

8000 

6500 

2.5 

2 

13000 

9000 

9 




19 

6000 

4600 

2.5 

5 

13000 

8000 

9 




30 

5000 

6000 

1.45 

3.4 

13000 

7000 

11 




47 

5000 

4400 

0.5 

1 

12000 

7000 

9 


1- The events are corrected for the galactic reddening effect (Schelegel et al. 1998). Additional host galaxy 
extinction is corrected for only in two objects, namely SNe 19901 and 1999em. 

2- The reported data corresponds to the low He I velocity case (see Sect. 3.2) 






























